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LIST OF SYMBOLS 


A, B,C,D = inertial parameters (slug-ft 2 ). 

C ai , C x ^, C X3 = hinge-couple components about axes parallel to the xy x 2 -, and x 3 -axes, respec- 
tively (ft-lb). 

C , Cy 2 ’ ( ~'y 3 = hinge-couple components about axes parallel to the clamshell body-fixed yy y 2 -, 
and y 3 -axes, respectively (ft-lb). 

Ci//, Cd, C<(> = cosines of the Euler angles i p, 6, and <j>, respectively. 

d 1 = clamshell hinge-axis displacement from the Xj-axis (the rocket vehicle longitudinal 
axis) (ft). 

d 2 = clamshell center of mass (c.m.) displacement from the x 1 x 2 -plane (the system bi- 
section plane) before clamshell deployment (ft). 

d 2 = clamshell c.m. displacement from the system base plane (ft). 

d 5 = clamshell c.m. displacement from the hinge axis (ft). 

d 6 = clamshell c.m. displacement from the x 2 * 3 -plane (the rocket vehicle system trans- 
verse plane containing its bary center) (ft). 

d n = clamshell c.m. displacement from the x^axis during deployment (ft). 
d 2 f = terminal value of d 2 (ft). 

F Xi , F X2 , F Xs = hinge-force components directed along axes parallel to the xy x 2 -, and x 3 -axes, 
respectively (lb). 

F , F y3 , F y3 = hinge-force components directed along axes parallel to the clamshell body-fixed 
y y y 2 -, and y 3 -axes, respectively (lb). 

F v F 2 , F 3 , F 4 , F 5 = inertial forces (lb). 


J c yi = \ m ( y\ + y\)dm 

J cy 2 = \ m ( y\ + y l )dm 


significant elements of the clamshell inertia matrix defined in 
terms of the clamshell body-fixed y-frame (slug-ft 2 ). 
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J cy 3 = ( iy\ + y\)dm I 

m I significant elements of the clamshell inertia matrix defined in 

| terms of the clamshell body-fixed y-frame (slug-ft 2 ) 

J c y rS m y ^ dm J 

J = rocket vehicle (minus clamshells) spin moment of inertia (moment of inertia 
about the Xj-axis) (slug-ft 2 ). 

J z i> J z ^, J Z3 = clamshell moments of inertia about the z t -, z 2 -, and z 3 -axes, respectively (slug-ft 2 ). 

K = direction cosine matrix. 

M yi , M y j M y3 = moments about the clamshell body-fixed y 2 -, y 2 -, and y 3 -axes, respectively (ft-lb). 

M Zi , M Z2 , M Z3 = moments about the clamshell body-fixed Zj-, z 2 -, and z 3 -axes, respectively (ft-lb). 

m = clamshell mass (slugs). 

p = position vector from the origin of the x-frame to the clamshell c.m. (ft). 

R l = component of the position vector from an inertial frame origin to the x-frame 
origin directed along the rocket vehicle body-fixed x^axis (ft). 

S\j/, SO, S<t> = sines of the Euler angles 0, and <p, respectively. 

t = elapsed time (s). 

jy = time at the end of the clamshell deployment phase and the beginning of the free- 
flight phase (s). 

U, V = momental parameters (ft-lb). 

JVj, W 2 , W 3 = clamshell free-flight rotational rate components about the Zj-, z 2 -, and z 3 -axes, 
respectively (s -1 ). 

= displacement vector for the ;'th point on the clamshell defined in terms of the 
X-frame (ft). 

j X £m ) = clamshell c.m. displacement vector defined in terms of the inertial X-frame (ft). 

{z ; J = displacement vector for the /th point on the clamshell defined in terms of the 
x-frame (the clamshell centroidal principal axis frame) (ft). 

a = angle between the XjX 2 -plane and the plane defined by the Xj-axis and the posi- 
tion vector p. 

(3 = angle in the clamshell mass-symmetry plane between the clamshell body-fixed 
y-frame and the clamshell centroidal principal axis set. 
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I 

I 

7 = clamshell roll-out angle (the angle between the x 2 -plane and theyjj^-plane). 

77 = angle between the XjX 2 -plane and the plane containing both the clamshell hinge 
axis and clamshell c.m. 

?7 0 = initial value of 77. 

ctp 7 p rif = terminal values of a, 7, and 77, respectively. 

\li, 6 ,<t> = Euler angles (see Figure 6 ). 

= rocket vehicle spin (rotational rate of the x-frame) (s -1 ). 
fly = terminal value of (s _1 ). 

J2 10 = initial value of (s -1 ). 
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INVESTIGATION OF A CLAMSHELL 
ROLL-OUT EJECTION CONCEPT 


by 

Lawrence F. Hatakeyama 
Goddard Space Flight Center 


INTRODUCTION 

In this report, a roll-out ejection concept for the release of clamshells from spinning sounding 
rockets is developed and discussed. The primary aim is to establish the desirability of this particular 
concept for use with sounding rockets. It will be seen that the conditions under which the clamshells 
are ejected impose requirements not considered in other applications. These requirements affect 
chiefly the manner in which the clamshells are to be ejected. 

At the present time, there are two general categories of ejectable payload-protection devices. The 
conceptually older and structurally simpler of these devices is the one piece nose cone. A nose cone is 
essentially a shell of revolution with its aft end faired and attached to the top stage of the rocket ve- 
hicle. The cone is tapered to a closed fore end. The payload is situated in the space bounded by the 
nose cone and the rocket vehicle. The nose cone is impelled at its ejection by springs or other means 
in the direction in which the rocket vehicle is pointed. Obviously, this is not attempted while the 
rocket vehicle is thrusting— it must occur under coasting conditions. If the rocket vehicle has a control 
system, it may be maneuvered so that the ejected nose cone does not present a collision hazard during 
a subsequent thrust phase. Unfortunately, sounding rockets do not now have this maneuvering capa- 
bility. Hence, the nose cone cannot be ejected safely until the sounding rocket is in its final coast 
phase and at an altitude where post ejection collision is not likely to occur. Thus, the performance of 
a sounding rocket can be impaired by its acceleration of excess mass. In the case of the Javelin (i.e., 
Honest John-Nike-Nike-X248) rocket vehicle, nose cone ejection is timed to occur 120 seconds after 
liftoff, when the vehicle is at an altitude of about 700,000 feet. At this point in the flight, any resid- 
ual X248 thrust and the drag deceleration difference between the ejected nose cone and rocket vehicle 
are considered to be negligible. It will be seen that this particular nose cone ejection is set for a time 
which occurs significantly later than is possible with clamshells. It can be seen also that, in general, 
nose cone ejection is troublesome when the payload is long and impossible when the payload compart- 
ment is bulbous. Guides or bumpers running the length of the payload and ejection actuators with 
long strokes are required in the former instance to avoid nose cone hang up. The nose cone can still be 
a source of trouble after it has cleared the rocket vehicle. It effectively continues to precede the 
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payload in the trajectory and may affect instrument readings by emitting particles and disturbing the 
environment in other ways. 

Clamshell systems are like nose cones, largely in overall shape. Each clamshell may be considered 
to be a longitudinal section of a shell of revolution. The clamshells are held together by bands, clamps, 
and the like, or they are attached to skin sections which can be ruptured at ejection time. On ejection, 
each clamshell is projected away from the longitudinal axis of the rocket vehicle; that is, its movement 
characteristically has a component which soon carries the clamshell out of the path of the payload. 
Therefore, the ejected clamshells do not continue to be collision hazards after they have cleared the 
rocket vehicle and payload. Thus, the time at which clamshell ejection is set to occur may be made 
meaningful in that it is not necessary to wait until the drag has dropped practically to zero before ejec- 
tion. In fact, a slight amount of drag will help to increase the longitudinal separation between the pay- 
load and ejected clamshells. Obviously, the application of rocket vehicle thrust can produce even 
greater separation. 

Ejection can occur for the Javelin when it is at an altitude of about 300,000 feet. At this point 
in its flight, it is about halfway into its X248 thrust phase. Perturbations due to X248 ignition and 
separation have been damped out, and the dynamic pressure has dropped to negligibly low values 
despite the considerable increase in vehicle velocity. The shape of the dynamic pressure profile for the 
Javelin is exemplified by the curve in Figure 1 . The X248 thrust phase occurs between the tick marks 
located at 56 and 97.9 s. In addition to permitting the recording of scientific data at lower altitudes 
than it was possible previously, clamshell ejection even at this point in the final boost phase will have a 
significant effect on the performance of the Javelin. Because about two thirds of the vehicle velocity 
at final burnout is due to the X248 thrust phase, the release of clamshells earlier in the flight can im- 
prove the vehicle’s performance (see Figure 2). 

Unfortunately, the conditions under which sounding rocket clamshells must operate have not 
been sufficiently considered in a number of designs. This situation may be partly due to the estab- 
lished success in the release of clamshells from nonspinning rocket vehicles, in which clamshells are 
disengaged and simply pitched out. It should be noted that in this case, the angular motion of each 
clamshell is restricted to rotation about that one of its centroidal principal axes normal to its mass 
symmetry plane. Hence, the motion of the clamshells during and after their deployment remains un- 
coupled and simple. This is not the case with clamshells pitched out from a spinning rocket vehicle. 
Instead, such an action causes each ejecting clamshell to rotate about all of its centroidal principal 
axes. The resulting complication greatly increases the extraneous tendencies of these clamshells and 
makes it difficult to design optimal constraints to control the clamshell motion during the deployment 
phase of ejection. 

The extraneous tendencies should be reduced, if not eliminated, by essentially limiting clamshell 
angular motion to rotation about a single principal axis. This requires that each clamshell be rolled 
out since it is rolling to begin with. This can be done by pivoting it about an axis through either its 
leading or trailing edge in the clamshell system sectioning planes. Extraneous rotational tendencies 
may yet be induced in the clamshell by reaction to constraints utilized to develop the desired rolling 
motion. Thus, care must still be exercised in the design of roll-out clamshell ejection mechanisms. 
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Figure 1— Dynamic pressure profile for a Javelin I 
(quadrant elevation angle) and carrying a 120-lb 






The angular momentum of a roll-out clamshell pivoted about its leading edge is increased on its 
deployment. This is caused by the displacement of its center of mass (c.m.) from the rocket vehicle 
longitudinal axis and the increase in its angular rate, which is a summation of its roll-out rate and the 
vehicular spin. The rocket vehicle will be despun by such clamshell deployment unless means are 
adopted to preclude it. This can be a troublesome endeavor since it tends to complicate the system 
design and increase its weight. 

A clamshell pivoted about its trailing edge, on the other hand, is despun on its deployment since 
its pivotal, i.e., roll-out, rate is in opposition to the rocket vehicle spin. Thus, the energy for clamshell 
deployment can be expected to come initially from the rotational energy of the clamshell. The rocket 
vehicle will experience a measure of despinning after the clamshell has pivoted to a given roll-out angle. 
This will definitely be the case when the clamshell is totally despun, i.e., when the magnitude of the 
clamshell roll-out rate equals that of the rocket vehicle spin. The clamshell may be disengaged at this 
point in its deployment to give its free-flight motion a purely translatory character. However, in an 
actual flight, it may be preferable to release the clamshell earlier, i.e., at a smaller roll-out angle to re- 
duce the extraneous torquing of the rocket vehicle during clamshell ejection. The selection of an opti- 
mum release angle is not obvious, particularly when too small an angle can result in a collision between 
the clamshells and the payload (this is the case with release at zero roll-out angle, i.e., instantaneous 
clamshell release). The effects of various system parameters must be investigated before any determi- 
nation can be made with respect to this aspect or any other aspect of this problem. 

ANALYTICAL ASSUMPTIONS 

The following analysis is concerned with the equations for the motion, forces, and couples gener- 
ated by clamshells released from spinning sounding rockets in accordance with a roll-out ejection con- 
cept. This concept requires that each ejecting clamshell be pivoted about an axis at its trailing edge in 
the system bisection plane so that its pivotal, i.e., roll-out, rate is in opposition to the rocket vehicle 
spin. Figure 3 illustrates the ejection sequence scheme viewed head-on to a rocket vehicle with a right- 
hand spin. 

In order to facilitate resolution of the problem, it is assumed that there is a problem symmetry 
which permits the characterization of the system dynamics by those of a single clamshell. Thus, the 



Figure 3— Trailing edge pivot type of roll-out clamshell system with right-hand vehicle spin. 
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clamshells are assumed to be dynamically matched, rigid bodies attached to a spinning rocket which is 
not coning in a significant manner when clamshell ejection is initiated. For convenience, it is assumed 
also that the damping and dissipative forces are negligible in comparison to the inertial forces. 

Three coordinate frames are used in the analysis of the deployment phase dynamics. One of 
these is the x-frame, which is centered at the vehicle system barycenter and orientated so that its 
Xj-axis is coincident with the rocket vehicle longitudinal axis and its x 2 -axis is directed in such a way 
that the clamshell system bisection plane is in the XjXj-plane. The x-frame may be assumed to be a 
rocket vehicle body-fixed frame since the barycenter may be considered to be stationary during the 
time required by clamshell deployment. The clamshell body-fixed y -fra me is centered at the clamshell 
c.m. and orientated so that itsyy, y 2 -, and y 3 -axes parallel the Xj-, x 2 -, and x 3 -axes, respectively, of the 
x-frame before clamshell ejection. The clamshell is constrained during its deployment to maintain the 
parallelism between the x 3 and y x axes. The z-frame is the clamshell centroidal principal axis set. It is 
oriented so that its z 2 -axis is normal to the clamshell mass symmetry plane and coincident with the 
y 2 -axis of the y-frame. The clamshell may be affixed with weights to rotate the Zj- and z 3 -axes and 
bring them into alignment with they^ andy 3 -axes, respectively, without changing the relationship be- 
tween the z 2 - and jyaxes. When this is done, the z-frame and y-frame are identical. Figure 4 illustrates 
the spatial relationship between the x-, y-, and z-frames. 

EQUATIONS FOR THE SYSTEM ANGULAR MOTION 

The equations for the angular motions of the system, derived 
by an application of Lagrange’s equation, may be written in a form 
suitable for digital computer solution as follows: 

& DUZBL 

1 AD-BC 

and 

.. _ AV-CU 
77 AD-BC ’ 

where 

B = - 2 (/ + rnd\) + 2 md $ d l cos rj , 

C-B, 

D = 2{J cyi + md\), 

U = -2(2mSl l rid s - mr] 2 d 5 )d l sin 77 , 

V = 2mnjd s d l sin rj , 

and 

AD - BC = 2J vx y cyi + md\) + 4 md 2 [J cyi + md 2 ( 1 - cos 2 17 )] 
^ 0 . 
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EQUATIONS FOR THE HINGE FORCES 

The hinge-force components are obtained by the application of Newton’s Second Law to the 
acceleration of the clamshell c.m. This yields 

^, = ^1 ’ 

F x = Fj suit? + (F 2 -F 3 ) cost? -F 4 sina -F 5 cosa , 
and 

F =Fj cost? - (F 2 ~F 3 ) sin T 7 +F 4 cos a-F 5 sin«, 

where 

Fj = mnd 5 , 

F 2 = mr) 2 d 5 , 


F 3 = 2mSl l r)d s , 

F 4 = m& x d n , 
and 

F s = mSlfa . 

EQUATIONS FOR THE HINGE COUPLES 


From Figures 4 and 5, it can be shown that 


and 


where 


and 


2 -Vi> 
C = M V + F d, + F d, , 

^3 ^3 J'l 1 ^2 3 

= 7W Zi cos (3 - A / Z3 sin (3 , 

M =M 7 , 
y 2 z 2 

M y3 = M Z3 cos j3 + M Zi sin (3 . 


c *2 = c y2 COS7 + c y 3 sin7, 

C x 2 = C y 2 cos 7 — C y 2 sin 7 ’ 

= F *2 COS7 - F X 3 Sin ^’ 

F y 3 = F x 3 cosy + F * 2 sin ^> 
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Solution of the preceding equations requires the 
application of Euler’s equation of motion to the 
problem; thus, 


^ = cos/3 , 

M. = (/ -/ )(J2j - 1 j) 2 cos/3 sin/3 , 

and 

M z 3 = ~ J z 3 ^i~^ sin 0 ’ 

where 

J, = J cos 2 /3 + / sin 2 /3 

z i cyi r cy 3 ^ 

- 2/ CJ , s cos/3 sin/3 , 

/ =/ , 

. z 2 ’ 
and 

/, = / sin 2 /3 + / cos 2 /3 

z 3 CJ*! ^ C^3 r 

+ 2 J cy5 cos/3 sin/3 . 

EQUATIONS FOR FREE FLIGHT 


Figure 5-Deployment phase of ejection showing one The f re e-flight displacements of the /th point 

ejecting clamshell. on c j ams h e n ma y b e expressed in terms of the 

X-frame, an inertial frame which is orientated so 

that its X t -, X 2 -, and X 3 -axes parallel the Xj-, x 2 -, and x 3 -axes of the x-frame at the instant of clamshell 
disengagement. This inertial frame translates at the rate established by the rocket vehicle at this time; 
thus, 


ra - m + k j 


where 


K = 


and 


cec<t> 

-se 

C6S<t> 

C\jj S6 C<j> + Sxjj S<p 

C\p C6 

C\l/S6S<j>-S\ljC<t> 

Sxfy S6 C<p - C\p S<p 

S^C8 

S\p S6 S(f> + Ci// C(j> 


^Ws sin r) f - & lf d lf sin ay >+ ld lf cos ay> . 


rjfd 5 cosTjj- + Slydy cos a. 


\dy sin ay. 


The square AT-matrix is a direction cosine matrix based on the Euler angle system shown in Figure 6. 
This angular system is a variant of a system used widely by aeronautical engineers. It is utilized to 
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2 



z 3 


Figure 6— Euler angles and free flight coordinate systems. 

simplify the determination of the initial Euler angles. From the construction in Figure 6, it can be 
shown that 


W 3 sin 0 + cos 0 
^ cos 0 ’ 

0 = cos 0 - IVj sin 0 , 

0 = W 2 + 0 sin 0 , 

0 =/ f , 
d = f‘ ddt , 

and 

0 = fidt-p. 

The z-frame components of the clamshell rotational rate may be obtained from Euler’s equations of 
motion for the free flight; thus, 
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and 


"i- T 


»WW 


W 2 = J 


z 2 




^ 3 = j 


z 3 


Wi =J W x dt + (fly-rif) cos 0 , 


W 2 = f W 2 dt , 

w 3 = J w 3 dt - (n if -T) f ) sin/ 3 . 


DISCUSSION 

Figures 7 through 25 illustrate the results of a study of a roll-out clamshell system for the Javelin 
rocket vehicle. The digital computer program and a data deck utilized in this study are listed in Ap- 
pendix A. Nominal system data, if the use of unaligned clamshells for which the |3-angle is not zero is 
assumed, are estimated to be as follows: 

J VXi = 7.5 slug-ft 2 , 

m = 0.4 slug , 

J cyi = 0.1 178 slug-ft 2 , 

/ = 0.3466 slug-ft 2 , 

/ = 0.4200 slug-ft 2 , 

J = 0.03219 slug-ft 2 , 
s 

d y = 0.8042 ft 2 , 
d 2 = 0.4286 ft 2 , 
d 3 = 1.456 ft 2 , 
flio = 9.5 rev/s , 

and 

^ = 515 ft-s' 2 . 
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When the study is applied to aligned clamshells, i.e., clamshells wherein the 0-angle has been zeroed, 
the applicable clamshell parameters are changed as follows: 

m - 0.4592 slug , 

/ = 0.1656 slug-ft 2 , 

J cy2 = 0.4654 slug-ft 2 , 

J = 0.5676 slug-ft 2 , 

c y 3 

J cys ~ 0.0 SlUg-ft 2 , 

d, = 0.3733 ft , 
and 

d 3 = 1.268 ft. 

These changes reflect the effects of alignment brought about by the attachment of two weights to each 
clamshell in a manner which results in minimum clamshell mass increase. 

The system motion and the hinge forces and couples generated by clamshell deployment are 
shown in Figures 7 through 1 0. There appears to be no significant difference between systems using 
unaligned and aligned clamshells according to these figures. 

It should be noted that the rocket vehicle is subject to slight spin-up followed by negligible de- 
spinning as the clamshells deploy. The individual and the total effects are of the order of a percent of 
the initial vehicular spin over the range of roll-out angles considered. No violation of angular momen- 
tum conservation is represented by the rocket vehicle spin-up because each clamshell is being despun 
as it rolls out. The vehicular spin-up signifies that the energy taken from the rotation of the clamshells 
is more than sufficient for their deployment. The excess energy is not large, so the spin-up is not sig- 
nificant. This observation applies also to the energy deficit which results in the rocket vehicle despin- 
ning at the larger roll-out angles. Thus, no special rocket vehicle despin avoidance devices are needed 
for the clamshell system simulated. 

Reversing the rocket vehicle spin permitted a comparative study of a roll-out system with pivot 
axis at the clamshell leading edge. As expected, such a system subjects the rocket vehicle to greater 
despinning and generates hinge forces and couples of considerably larger magnitudes than the system 
with trailing edge pivot. These effects, illustrated in Figures 1 1 and 12, are attributed to the fact that 
the clamshells are spun up as they roll out. It will be seen that this spin-up also raises the minimum 
roll-out angle at which the clamshells can be safely disengaged. Thus, a roll-out system with pivot axis 
located at the clamshell trailing edge is preferable to a system with pivot axis at the leading edge. 

It may be inferred from Figures 8 and 1 0 that the hinge couples are more significant to the sys- 
tem designer than the hinge forces. Thus, Figures 1 3 through 1 6 are included to illustrate the effects 
of rocket vehicle spin and longitudinal acceleration on C and C , the hinge couples which oppose 
the clamshell pitching and yawing tendencies, respectively. As expected, the rocket vehicle spin at the 
higher levels investigated produces a decidely bad effect on C x and C Xj . On the other hand, the rocket 
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vehicle longitudinal acceleration tends to reduce the maximum magnitude of C by shifting its time 
trace upward. No such beneficial effect is incurred for C despite a similar upward shifting of its time 
trace. Whatever the case may be, the magnitudes of C x and C indicate that serious consideration 
should be given to reducing the rocket vehicle spin to al»out a half of that presently utilized. Use of a 
lower rocket vehicle spin can improve the X248 motor performance in addition to moderating the de- 
sign requirements of the clamshell system. 

The X 2 X 3 projections of the near free-flight displacements of clamshells released at roll-out angles 
of 12.5, 15, 30, and 60 deg are shown in Figures 17 through 23. Except in Figure 19, these projec- 
tions are for aligned clamshells. Since the X 1 X 2 and X t X 3 projections for these clamshells are straight 
lines, and therefore of little interest, they are not presented. Figure 19 shows that the near free-flight 
displacements of an unaligned clamshell under the conditions considered is not markedly different 
from that of an aligned clamshell. It is possible that conditions beyond the scope of this study could 
produce effects requiring further investigation. 

Figures 17 and 20 show that roll-out clamshells can be released too soon. In each case, the clam- 
shell rotational magnitude is too high for release at the roll-out angle shown. Obviously, clamshell re- 
lease can take place safely at a lower roll-out angle with the trailing edge pivot type of system because 
the clamshells are subject to despinning and the offending parts are displaced farther away from the 
payload when disengagement occurs. The rotation of each clamshell will be near zero, and its free- 
flight motion thereby will be almost purely translatory when the clamshell is released at a roll-out 
angle of 60 deg. This effect occurs near 60 deg for the system under consideration at the various 
vehicle spin rates shown in Figure 24. Indeed, the angular motion of the system can be characterized 
by the reduced forms contained in Figure 25. This figure shows that the relationship between and 
rj is constant for any given roll-out angle. The locus of points in Af-space through which a given part of 
the clamshell passes is fixed therefore by the 7 -angle at which disengagement occurs. The vehicular 
spin merely affects the rate at which such a given set of points in A-space is traversed. Thus, clam- 
shells which can be released safely at 1 5 deg when vehicular spin is 9.5 rev/s can also be released safely 
at this angle at any other positive vehicular spin if the system can bear the loads imposed upon it. That 
is, clamshell release for the system under consideration can be programmed for roll-out angles between 
1 5 and 60 deg. Choice of the lower angles will be influenced by the desire to reduce unbalanced torqu- 
ing of the rocket vehicle during clamshell deployment. This torquing may arise from vehicular coning 
motion, clamshell mismatch, and the “yo-effect” caused by nonsimultaneous release of clamshells. On 
the other hand, release at a higher roll-out angle is desirable because it results in the ejection of clam- 
shells with reduced rotational motion and lowered likelihood of collision with the payload. 

The system angular motions are not affected by the rocket vehicle longitudinal acceleration. The 
system characteristics discussed in the preceding paragraph will therefore be independent of deviations 
in rocket vehicle thrust. Since the various acceleration levels are normally associated with different 
system mass properties, it was expected that the curves in the figures discussed would reflect this fact. 
This mass effect, however, tends to be a minor one since it involves the interchange of a relatively 
small amount of energy between the rocket vehicle and the deploying clamshells. 
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Figure 9— System motion with aligned clamshells. 


Figure 10— Forces and couples with aligned clamshells. 
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Figure 19-Near free-flight displacements of an unaligned roll-out clam- 
shell released at 15 deg. 



10 12 

X 3 (ft) 


Figure 20— Effect of rocket-vehicle spin reversal on the near free-flight 
of an aligned clamshell released at 15 deg. 
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Figure 21— Near tree-flight displacements of an aligned roll-out clamshell 
released at 30 deg. 



Figure 22-Near free-flight displacements of an aligned roll-out clam- 
shell released at 60 deg. 
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Figure 23— Effect of rocket-vehicle spin reversal on the near free-f light 
of an aligned clamshell released at 60 deg. 



Figure 24— Clamshell rotational rate. 
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Appendix A 

Source Listing of Program "ROC" and a Data Deck 
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n o o 


K> 

4 ^ 


// EXEC F0RTRANG,PARM=* NA ME= ROC , DECK 1 
//SOURCE. SYSPUNCH DD 0 SN=&DECK , SYSOUT=B 
//SOURCE. SYSlN DD * 

C 2/21/69 - L.F.H. 

MOD. 8/21/69 - L.F 


• H • 


ROCOOOIO 
ROC00020 
ROCOO03O 

...MAIN PROGRAM 'ROC 1 ... ROCOOOIO 

PROG. ROC AND ITS SUBSIDIARY SUBPROGRAMS MAY BE UTILISED TO STUDY ROCOOOSO 
ROTH THE DEPLOYMENT AND THE FREE FLIGHT PHASES OF CLAMSHELL EJECTION. RDC00060 
THE TYPE OF SYSTEM THAT CAN BE STUDIED IS BASED ON A UNIQUE CLAMSHELL ROC0OO7O 
ROLL-OUT EJECTION CCNCFPT. THE FREE FLIGHT PHASE OF A GIVEN CASE IN AROC00080 
JOB RUN MAY BE SKIPPED IF DESIRED. THE DATA LOADING IS SET SO THAT R0C00090 
THE TERMINAL CARD IN EACH LOGICAL SUBSET OF CARDS REPRESENTING A CASE ROCOOlOO 
INPUT 8E PUNCHED WITH AN INTEGER OF THE FORM ' LLLKJ 1 AND ADJUSTED AS ROCOOllO 
DESIRED BETWEEN COLUMNS l THROUGH 8, INCLUSIVE. IF f J' IS ZERO , THE R0C00L20 
FREE FLIGHT PHASE OF THE CASE IS SKIPPED. IF • K* IS ZERO, THE CASE ISR0C00130 
THE LAST IN THE JOB RUN TO BE PROCESSFD. THE CASE NO. IS OPTIONALLY R0C00140 
ENTERED BY PUNCHING UP TO THREE DIGITS IN • LLL 1 ; IF IT IS OMITTED, THERDC00150 


CASE NO. 


IS EQUAL TO THE PRECEDING CASE NO. PLUS ONE. 

L.F.H. 8/21/69. 


120 FORMAT ( 1 0 * , 39X,*.«. J1 = '.II, 1 ...' ) 

130 FORMAT C 80X,'RTIME = ',F8.3,* SEC . * ) 

EXTERNAL SFTl , AUXl , DERI , AUX2 ,OUT 1 , 

X SET 2, AUX3,DER2, AUX4, OUT 2 


COMMON 

REAL*8 


/ LFH2 / 


COMMON / DATA / 
EQUIVALENCE 


NX ( 47 ) ,CBLOCKl 20,220) 

Z (250) , 

H,X( 3) , Y( 8, 3) , 

DER ( 8 , 4 ) 

Z ,NEQl,NEQ2,Jl.J2,J3fJ4,J5 


1 

2 

300 CALL ST l ME 
CALL LOAD ( 
IF( NX ( 2 0 ) 
KK = NX ( 9 ) 


H, Z ( 44 ) ),( X ( 1 ) , Z ( 6 1 ) ),( Y ( 1 , 1 ) , Z ( 64 ) ), 
DER I 1 , l ) , Z ( 9 1 ) ) 


NX( 22) 
NX ( 44 ) 


. NE • 0 ) GO TO 990 


JJ = MOD ( KK , 1 0 ) 

JK = MOD ( KK, 100 )/l0 
JL = MOD ( KK, 100000 )/100 
I F ( JL .GT. 0 ) J4 = JL 

...1ST ( DEPLOYMENT ) PHASE. 

CALL NIT 1 

CALL RK ( X, Y,DER, 3,H, Jl, J2, 

X SET1, AUXl, DERI, AUX2,0UT1 ) 

I F ( Jl .GT. 4 ) GO TU 500 
CALL AUXl 
CALL DER1( 

CALL AUX 2 
CALL OUT 1 ( 


X , Y , DER , 3,1 ) 


1 ) 


RGC00160 
ROC00170 
RQC00180 
R0C001 90 
R0C00200 
ROC0021 0 
R0C00220 
ROC002 30 
R0C00240 
ROC00250 
RQC00260 
R0C00270 
R0C00280 
R0C00290 
R0C00300 
R0C00310 
ROC00320 
R0C00330 
R0C00340 
R0C00350 
R0C00360 
RQC00370 
R0C00380 
R0C00390 
R0C00400 
R0C004 10 
ROC00420 
R0C00430 
ROC00440 
R0C00450 
R0C00460 
ROC 00470 
R 0C00480 
ROC 00490 
ROC00500 
ROC 005 10 
RGC00520 



n o o o n o o o o ooo 


to 

Oi 


I F t JJ . EQ. 0 > GO TO 550 

..•2ND ( FREE FLIGHT ) PHASE... 


CALL N I T 2 

CALL RK( X,Y,DER,6,H,Jl,4, 

X SET2,AUX3,DER2,AUX4,0UT2 ) 

I F ( JL .GT. 4 ) GO TO 500 
CALL AUX 3 

CALL DER 2 ( X,Y,0ER,6,l ) 

CALL AUX4 
CALL OUT 2 { 1 ) 

GO TO 550 

500 WRITE! 6,120 ) Jl 
550 CALL T TI ME ( JTI ME ) 

RT I ME = FLU AT l JTIME*26 J/1000000. 
WRITE! 6,130 ) RTIMF 
IF! JK .NE. 0 ) GO TO 300 
GO TO 999 

990 IF! JK .EQ. 0 ) GO TO 999 
GO TO 300 

999 STUP 
END 

2/22/67 - L.F.H. 

MOD. 10/23/68 - L.F.H. 

SUBROUTINE LOAD 


...PROD. VERS... 


C 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


110 FORMAT! 20A4 ) 

120 FORMAT! 1 X, I 2, 5X, 2 044, 2 X, Z 8 , 1 X ,Z 8 , 1 X , l 8 ) 
130 FORMAT! 15X,'...RFAD ERROR IN LOAD... 1 ) 


INTFGER*2 
COMMON / LFH2 / 

1 

2 

3 

4 


K 3 , M 1 

FWD,RD,RZ(4). 

Rl,R2fR3,R4,R5,R6,R7,R8,R9,RY 13) , 
N 1 , N ? , N 3 , N4 , 

CARO ( 20 1 , BLANK , 

R30,PSW{ 7 ) , K3 , Ml , COLOCK (20,2201 


F WO IS NOT USED BY SUBP. LOAD OR SUBP. EDATA. 

RD CONTAINS SLOP. LOAD'S SAVE ARtA ADDR. 

RZ'S ARE NOT USED BY SLBP. LOAD CR SUBP. EDATA. 

R 1 CONTAINS THE RETURN ON AN SPIE ISSUED IN SUBP. EDATA. 

R2 CONTAINS A(OUST), THE SHIFTED B- ADDR • USFD IN SUBP. EDATA. 

R 3 CONTAINS THE A ( CURRENT LOAD POINT) UN EACH RETURN FROM SUBP. 
EDATA. 

R4 CONTAINS A ( LFH? ) , I.F., SUBP. EDATA'S SAVE AREA ADDR. 

R 5 CONTAINS A(CARD), T.E., A(LFH2<-88). 

R6 CONTAINS A { CURRENT CARD CHARACTER IMAGF). 

R7 CONTAINS A t CARD + 79 ) , T.F., MLFH2U67). 


R0C00530 
R0C00540 
R0C00550 
R0C00560 
ROC005 70 
ROC 00580 
ROC 00590 
ROC00600 
R0C006 10 
ROC00620 
RDC00630 
R0C00640 
ROC00650 
ROC00660 
R0C00670 
R0C00680 
ROC00690 
ROC00700 
ROC0071 0 
R0C00720 
R0C00730 
R0C00740 
R0C00750 
ROC00760 
ROC 00770 
R0C00780 
LOADOOIO 
L0AD0020 
L0AD0030 
L0AD0040 
L0AD0050 
LOAD0060 
L0AD0070 
L OAD0080 
L0AD0090 
LUAD0100 
LOADOLLO 
LOADO 120 
LO ADO 1 30 
LOADOI 40 
L0AD0150 
LQAD0L60 
LOADOI 70 
LOADO 180 
LOADO 190 
L0AD0200 
L0AD0210 
L0AD022 0 
L0AD0230 
L0AD0240 
L0ADO250 
L0AD0260 
LQAD0270 
L0AD0280 
LOAD029O 



ooo ooononooonoooonnooonnnon 




R8 CONTAINS A t LOB ) , AN ADQR. INTERNAL TO SUBP. EDAT A. 

R9 CONTAINS At 3) . 

RY • S ARE NOT USED BY SUBP. LOAD OR SUBP. EDATA. 

N1 IS A BRANCHING CONTROL SET BY SUBP. EDATA FOR SUBP. LOAD. 
N2 IS A CARD COLUMN FRROR INDICATOR SET BY SUBP. EDATA. 


L0AD0300 
LOAD0310 
L0AD0320 
LOAD0330 
LOADO 340 


N3 IS A COUNT OF CAROS •READ* BY SUBP. LOAD AND PROC. BY SJBP. ED AT ALG4D0350 

N4 = 32767 IS SET BY SUBP. LOAD TO SIGNAL A •READING* FRROR. L0AD0360 

CARD ( 2 0) CONTAINS THE CHARACTERS OF THE CURRENT CARD IMAGE 1 READ* BYLG AD0370 
SUBP. LUAD AND PROC. BY SUBP. EDATA. L0AD0380 

R30 CONTAINS AtINIT. LOAD POINT). LOAD0390 

PSW CONTAINS THE O-PSW WHENEVFR SUBP. EDATA IS SUBJECTED TO A PROG. L0AD0430 
INTERRUPTION WHICH IS CHANNELED BY THE * $P I E * IN SUBP. EDATA. LOA00410 
K3 = 1 IS SET BY SUBP. LOAD TO SIGNAL AN ENCOUNTER WITH AN EOF. L0AD0420 

Ml = 1 IS SET BY SUBP. LOAD AFTER ALL CARDS IN A FILE HAVE BEEN LDAD0430 

PROC. BY SUBP. EDATA. LOADCmO 

CBLOCK IS THE AREA WHERE A BLOCK UF CARD IMAGES ARE ’READ' IN BY L0AD0450 
SUBP. LOAD. THIS BLOCK MAY CONTAIN J1 CARO IMAGES OR ALL CARD L0AD0460 
IMAGES UP TO AN 'EOF*, WHICHEVER IS LESS. WHEN THE NO. OF L0AD0470 

CARDS IN A FILE IS LESS THAN Jt * SUBP. LOAO WILL 'READ* ALL OF L0A00480 
THE CARDS IN THF FILE AND ENCOUNTER THE *EOF*. IN SUCH AN L0AD0490 

EVENT, SUBP. LOAD WILL SET K3 = 1 AND Ml = 0. AFTER THE LAST L0A00500 
CARD IN THE FILF HAS BEEN PROC. BY SUBP. EDATA* SUBP. LOAD L0A0051G 

SETS Ml = 1. LGAD0520 

L0AD0530 

DATA R / 4HRR. • /, L0AD0540 

1 J1 / 220 /, LOAD0550 

2 K2 / 0 / L0AD0560 

LOAD0570 
L0AD058O 
LOAD0590 
L0AD0600 
LOADO 510 
LOAD0620 
LOAD0630 
L0AD0640 

IF( K2 .LT. 1 ) GO TO 200 LOAD0650 

Kl = K2 + 1 LOAD0660 

GO TO 320 L0AD0670 

LOAD0680 
LOA00690 
L0AD0700 
LOAD07 10 
L0AD0720 
LOAD0730 
L0ADD740 
LOAD 0750 
L0AD0760 
L0AD0770 

) GO TO 470 L 0 A DO 7 8 0 

L0AD0790 
L0AD0800 
LOADO 810 
L QAD0820 
L0AD08 30 
LOADD8<tO 


NOTE. • .CBL0CK(20»M) REQUIRES THAT J1 .EQ. M. 

Nl = 1 
N2 = 0 
N3 = 0 
N4 = 0 


200 DO 300 J = l,Ji 
300 CBLOCK 1 1 , J ) = R 
Kl = 1 
K3 = 0 
Ml = 0 

READt 5* 1 10 , END=600 , ERR= 700 ) CBLOCK 


320 DO 400 K = K 1 , J 1 

I F ( CBLOCK ( 1 * K ) .EQ. R 
K2 = K 

330 DO 350 J = 1*20 
350 CARD(J) = CBLOCK ( J , K ) 
N3 = N3 + 1 



on n no non 


CALL EOATA 

GO TO { 400,450,370 ) ,Nl 
C 

370 I F ( N2 • E0. 0 ) GO TO 400 

I F ( PSW(l) .NE. 0.0 ) PRINT 120, N2 , CARD, PSW, R4 
I F ( PSW(l) .EC. 0.0 ) PRINT 120, M2* CARD 
400 CONTINUE 

I F ( K3 .LT. 1 ) GO TO 200 
C 

450 I F ( K2 .LT. J1 ) GO TO 500 
470 K2 = 0 

1 F l K3 .EG. 1 ) Ml = 1 
500 RETURN 
C 

600 10=1 

GO TO 320 
C 

700 N4 = 32767 
PRINT 130 
GO TO 500 
END 

C 2/18/69 - L.F.H. 

C MOO. 7/07/69 - L.F.H. 

C 

SUBROUTINE NIT1 

...DATA INIT. FOR THE 1ST PHASE... 


RE AL*8 

1 

2 

3 

4 

5 

COMMON / OATA / 
EQUIVALENCE 

1 

2 

3 

4 

5 


Z1250) , 

T 0 ,H0 , Q ( 8 ) , 

DU0),JCY(6),JZ(3), 

H, ETAZ, BETA, X,Y(8), 

DC ( 5 ) , DS ( 5 ) , 

TWPI ,CRTD,CDTR,P1D2, PID4 
Z,NE0l,NE02,JlfJ2,J3,J4,J5 

( TO , 7 ( 1 ) ),( HO , Z { 2 ) ),( Q ( 1 J , Z ( 3 ) ), 

( 0 ( 1 ) , Z ( 1 1 ) ),( JCYIl), Z(?2) ),( JZ ( 1 ) , Z ( 28 ) 
( H , Z ( 44 ) ),( ETAZ,Z(46) ),( BETA,Z(47) ), 

{ X , Z ( 6 1 ) ),( Y C 1 ) , Z ( 64 ) ), 

< ocm,z(i3n ),( osm,zu36) ) 


COMMON / CONS / T WP I , CRTO , CDTR , P I D? , P I D4 


X = TO 
H = HO 

ETAZ = DAT AN ( 0(2)/0{l) > 

0(3) = E TAZ *CR TO 
Y ( 1 ) = Q ( 1 ) *TWP I 
Y ( 2 ) = Q(2)*TWPI 
Y ( 3 ) = ETAZ 

D ( 5 ) = DSQRT ( D ( 1 ) **2 + D(2)**2 ) 


...SET UP BETA AND THE PRIN. M. OF INERTIAS... 


L0AD0650 
LOA00860 
L0AD0870 
LOAD0830 
L0AD0890 
LOAD0900 
LOAD0910 
L0AD0920 
L0AD0930 
L0AD0940 
L0AD0950 
LQAD0960 
LOAD0970 
LOAD0980 
LOAD0990 
LOADIOOO 
LOAD 1010 
LUAD1020 
L0AD1030 
L0AD1040 
LOAD! 050 
1NIT0010 
INI T 0020 
1NIT0030 
I N I TO 040 
IN1T0050 
INIT0060 
1NIT0070 
INIT0080 
1NIT0090 
1NIT0100 
INIT0110 
1NIT0120 
INITO 130 
INI T0140 
1 NIT01 50 
1NIT0160 
) , LN I TO 170 
1NIT0180 
1NIT0190 
INI T0200 
1NIT0210 
1NIT0220 
INIT0230 
INIT0240 
INIT0250 
INI T0260 
1NIT0270 
1NIT0280 
INIT0290 
1N1T0300 
1 N I TO 310 
IN l TO 3 20 
1NIT0330 
1NIT0340 



Is) 

00 


IF( JCY(l) • EQ. JCYO) ) GO TO 300 
BETA = 0 • 5D0+DAT AM ( 2.000* JC Y( 5) /( JCY(3) 
I F ( BETA .GT. PID4 ) BETA = BETA - PID2 


JCY(l) ) I 


I F { BETA + PID4 .LT. 
GO TO 310 
300 BETA = O.ODO 

DCOS ( BETA ) 


O.ODO ) BETA 


BETA + P l 02 


310 Dcm 
DS ( A ) 
JZ( l) 

1 

J Z ( 2 ) 
JZ ( 3 ) 
l 


DS I N ( BETA ) 

JCY(1)*DC(4)**2 + JCY(3)*0S(4)**2 

- 2.0D0*JCY(5)*0C(4)*DS14) 

JCY ( 2 ) 

JCY ( 1 ) * DS ( 4 ) **2 + JCY ( 3 ) *DC ( 4 ) **2 

+■ 2.0D0*JCY( 5)*QC(4)*DS(4) 


C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


990 RETURN 
END 

2/14/69 - L.F.H. 

MOO. 8/22/69 - L.F.H. 

SUBROUTINE RKl / X/ , /Y/ , /DER/ , /NEC/ , /H/ , / J 1 / , / JF/, 
l SETH, AUX1,DERI V, AUX2, OUT ) 

. . .RUNGE-KUTTA 4TH ORDER I NTEGR A TOR. . . 

FIXEO STEP INTEGRATION EXCEPT FOR THE TERMINAL PROCEDURE. 

X - INDEPENDENT VARIABLE. 

Y 1 S - DEPENDENT VARIABLES. 

DER'S - DERIVATIVES OF THE Y*S. 

NEC) - NO. OF DERIVATIVE EQUATIONS. 

H - INTEGRATION STEP SIZE. 

Jl - BRANCHING PARAMETER SET BY SUBP. SETH FDR SUBP. RK. 

JF - OUTPUT PRINT FREQUENCY. 

SUBP. SETH - ADJUSTS H DURING THE TERMINAL INTEGRATION PROCESS AND 
SETS Jl = 1,2, 3, 4, OR 5 DEPENDING ON WHETHER H IS LEFT UNCHANGED OR 
ADJUSTED WITHOUT THE NEED TO RESET X(l) AND Y(K,1>, THE INTEGRATION 
PROCESS IS TO BE ENDED, H IS CHANGED AND THE PRECEDING INTEGRATION 
STEP IS TO BE REPEATED WITH THE NEWER H-VALUE AND THE RESETTED X AND 
Y-VALUGS, H MUST BE CHANGED BEYOND THE MAX. ALLOWED NO. OF TIMES, OR 
AN ABNORMAL END OF THE INTEGRATION PROCESS IS RFQU I RED ,RE SPE C TI VE LY. 
SUBP. SETH DEPENDS ON J1=0 WHENEVER SUBP. RK IS CALLED. 

SUBP. AUXl - COMPUTES, l.E. UPDATES, THE RFQ'D OATA FOR SU3P. DERIV 
SUBP. DERIV - COMPUTES THE DERIV 1 S OF THE Y* S AT THE SUBSTEP POINTS 
SUBP. AUX2 - COMPUTES, l.E. UPOATES, THE ADDITIONAL DATA FOR OUTPUT 
SUBP. OUT - IS THF OUTPUT SUBPROGRAM. 


REAL+8 


Jl = 0 
JN = 0 


XI 3),Y(NE0,3) ,DFR(NEQ,4) ,H, 
CHH,HD6 


400 X ( 2 ) = X ( 1 ) 

DO 410 K = 1 , NEQ 


1NIT0350 

INI T0360 

1NIT0370 

1 N I TO 380 

INI T0390 

IN I T 0400 

1NIT0410 

INI T0420 

1N1T0430 

INI T0440 

INI T0450 

INI T0460 

1NITQ470 

1NIT0480 

1NIT0490 

INIT0500 

RKOOOOIO 

RK000020 

RK000030 

RK000040 

RK000050 

RK000060 

RKOOOO 70 

RK000080 

RK000090 

RKOOOIOO 

RK00011 0 

RK000120 

RK 0001 30 

RK000140 

RKOOO 150 

RKOOOl 60 

RK000170 

RK000180 

RKOOOl 90 

RK000200 

RK0002 10 

RK 000220 

RK000230 

RK000240 

RK000250 

• RK0002 50 
. RK000270 

• RK 000280 
RK000290 
RK000300 
RK000310 
RK000320 
RKOOO 330 
RK000340 
RK000350 
RKOOO 360 
RK000370 
RKOOO 380 
RKOOO 3 90 


o o o n 


to 

VO 


410 Y(K,2) = Y(K,l> 

IF{ J1 .NE. 3 ) GO TC 415 
I F { JN .GT. 1 -AND. -Moot JN-l,JF ) 

C 

415 CALL SETHI X , Y , DER ,NEC f H , J L ) 

GO TO ( 420,990,400,500,500 >,Ji 
420 CONTINUE 
C 

CALL AUX1 

CALL DERIVt X , Y , OER , NEC, 1 ) 

IF! MODI JN, JF ) .NE. 0 ) GO TO 425 
CALL AUX 2 
CALL OUT ( JN ) 

425 CONTINUE 
C 

CHH = 0.5D0*H 
C 

DO 440 J = 2,4 
I F ( J .EQ. 4 ) CHH = H 
X(l) = X I 2 ) + CHH 
DO 430 K = l,NFQ 

430 Y ( K , 1 ) = Y ( K ,2 ) + CHH+DER ( K, J-l ) 
CALL AUX 1 

CALL DERIVC X , Y, DER , NEO, J ) 

440 CONTINUE 


RK000400 

RK000410 

NF. 0 ) JN = JN - 1 RK000420 

RK000430 
RK000440 
RK000450 
RK000460 
RK000470 
RK000480 
RK000490 
RK000500 
RK0005L0 
RK000520 
RK 000530 
RK000540 
RK000550 
RK 000 5 60 
RK000570 
RK000580 
RK000590 
RK000600 
RK00061 0 
RK000620 
RK000630 
RK000640 
RK000650 
RK000660 
RK000670 
RK000680 
RK000690 
RK 000700 
RK0007 10 
RK000720 
RK000730 
RK000740 
RK000750 
RK000760 
RK000770 
RK000780 
RK000790 
1S6T0010 
1 SET0020 
1SET0030 
l SET 0040 
1SET0050 
i SE T0060 
ISET0070 
LSET0080 
LSET0090 
1SET0100 
lSETOHO 
1 SET0120 
i SET0130 
1SET0140 
ISET0150 


HD6 = H/6.0D0 
X ( 3 ) = X I 2 ) 

DO 470 K = l , NEQ 
Y ( K , 3 } = Y { K , 2 ) 

Y ( K , 1) = DER ( K , I ) + 2 • ODO* ( DER ( K ,2 ) 
470 Y ( K , 1 ) = Y { K , 2 ) + HD6*Y(K,L) 


+ DER ( K , 3 ) ) + DER ( K , 4 ) 


JN = JN + 
GO TO 400 


1 


500 CONTINUE 


990 RETURN 
END 

2/14/69 - L.F.H. 

MOD. 8/22/69 - L.F.H. 

SUBROUTINE SETH /X/ , /Y/ , /OER/ , /NEC/ , /H/ ,/ Jl / } 


...TERMINAL H CONTROLLER FUR THF 1ST PHASE... 

X(l) AND Y ( K , 1 ) ARE RESET WHEN AN l NT EG. STEP IS TO BF REPEATED. 


RFAL*8 


COMMON / DATA / 
EQUIVALENCE 


7 ( 250 ), 

YT,YRPS, 

FTA? V 0ER(MEQ,4) » 
H,X(NEQ),Y(NcQ, 3) , 
GAMD,TWPI,CRTD,CDTR,P ID2,PID4 
If KKK( 7) 



u> 

o 


X 

X 


c 

c 

c 

c 

c 


c 

c 

c 


c 


c 


c 


c 


c 

c 

c 

c 

c 

c 


COMMON 
I F { Jl 


( YT ,Z(42) ),( YFPS, 7.(43) 

I ETAZ,2(46) ),( J3 f KKK(5) 

/ CONS / TWPI ,CRTD',CDTR,PID2,PID4 

.EQ. 0 ) JC = 0 

• • .CHECK FOR TERMINAL CONDITIONS... 


) 


GAMD = i Y( 3,1) - ETAZ )*CRTD 
IF! GAMD .LT. O.ODO ) GO TO 500 

I F ( DABS ( ( GAMD - YT )/YT ) .LE. YFPS ) GO TO 200 
I F ( GAMD .G T. Y T ) GO TO 300 
I F ( J1 .EQ. 0 > GO TO 100 

I F ( Y ( 2 , 1 ) .LE. 0.000 -AND. DER(2,1) .LE. O.ODO ) GO TU 500 


...SET J1 AS RFQUI RED • • • 


100 Jl = 1 

GO TO 990 

200 Jl = 2 

GO TO 990 

300 JC = JC + 1 

I F ( JC .GT. J3 ) GO TO 400 
I F ( Ji .EQ. 0 ) GO TO 500 
H = 0.5D0*H 
X(l) = X ( 3 ) 

DO 330 K = 1 , NEQ 
330 Y(K,l> = Y ( K * 3 ) 

Jl = 3 
GO TO 990 

400 Jl = 4 

GO TO 990 

500 Jl = 5 
990 RETURN 
END 

2/14/69 - L.F.H. 

MOD. 7/07/69 - L.F.H. 

SUBROUTINE AUX1 

...DATA UPDATER FOR THE 1ST PHASE DERIV. SUBP... 

REAL*8 Z ( 250 ) t 

1 0(10) ,M, JCY<6) , 

2 JVXl,GAM,ETAZ?Y(8), 

3 DC ( 5 ) , DS ( 5) , 

4 F ( 6 ) * F23* 

5 QOM.QET, 

6 A , B, C , DO , U, V 

COMMON / DATA / Z ,NEQl , NEQ2 , Jl, J2, J3, J4, J5 


l S6 TO 1 60 
1SFT0170 
1SET0130 
1SET0190 
1SET0200 
1 SET 02 1 0 
ISFT0220 
1 SET0230 
1 SET 0 2A0 
1 SET 0250 
1SET0260 
ISET0270 
1 SET0280 
1SET0290 
ISET0300 
1SET0310 
1SET0320 
1 SE TO 330 
1 SET0340 
1SET0350 
ISET0360 
1 SET0370 
1 SET0380 
1 SET 0 390 
1 SET0400 
1SET0410 
1 SET 0420 
1 SET0430 
1SET0440 
1SFT0450 
1 SET 0460 
1 SET0470 
1SET0480 
1 SET 0490 
1 SET0500 
1SET0510 
1 SET0520 
1 SET0530 
1 SET0540 
1 SET0550 
1 AUX0010 
I AUX0020 
l AUX0030 
1AUX0040 
1AUX0050 
1 AUX0060 
1AUX0070 
1 AUX0080 
1 AUX0090 
lAUXOiOO 
1AUX0110 
1AUX0120 
1AUX0130 
1AUX0140 
1AUX0150 



non no o o o on o oo 


EQUIVALENCE 

1 

2 

3 

4 

5 

COMMON / D3XI 


( nm,z(in ) f < m,z( 2 d i,( jc v 1 1 ) , z <221 ), 

( JVXl ,7(31) ),( GA M, Z ( 32 ) >,( ETAZ f Z<46) ), 

( Y ( 1 ) * Z ( 64 ) ), 

( DC ( 1 ) » Z ( 1 3 1 ) ),( DS ( l ) » Z ( 1 36 ) ), 

( F(l) ,7(141) ), ( F23, Z ( 147 ) ) 

A, B,C, ODi U, V 


DC U ) = DCOS ( Y ( 3 ) ) 
D$U) = DS I N ( Y(3) ) 


GAM = Y C 3 ) - ETAZ 
OC ( 2 ) = DCUS ( GAM ) 

D$ ( 2 ) = DS! N ( GAM ) 

D ( 7 ) = DSQRTt 01i)**2 + Di5)**2 - 2 . 0D0*D l 5 ) *011 ) *DC U) ) 

F(2) = M* ( Y<2)**2 )*0(5) 

F ( 3 ) = 2 *ODO*M*Y (1)*Y(2)*D(5) 

F23 = F ( 2 ) - F t 3 ) 

CALL FF6 

DD = 2. 0D0*( JCYU) + M*{ D(5)**2 ) ) 

A = JVXl + 2 • ODO* ( JCYU) + M*{ D(7)**2 ) ) 

B = - DD + 2«0D0*M*DI5)*D11)*DCC 1) 

C = B 


OOM = O.ODO 
QErT = F I 6 ) *D ( 1 ) 

U = QOM + 2.0D0*F23*D(1)*DS(1) 

V = GET + 2 • ODO*M* ( Y<1)**2 ) *DC 5 ) *D t l ) *DS ( 1) 


990 RETURN 
END 

2/14/69 - L. F *H. 

SUBROUTINE 0ER1< X f Y , DFR , NEG, J ) 


• • • DER I V • FQ* $• FOR THE 1ST PHASE... 


C 

C 

c 


c 


REAL *8 
l 


COMMON / DBXl / 

DET = A*D - B*C 

DER (It J) = ( D*U 
DER ( 2 , J ) = ( A*V 
DER t 3, J ) = Y ( 2 ) 


990 RETURN 
END 

C 2/14/69 - L.F.H. 


X,Y(NEQ) ,DER{NEQ,4) , 
A,B,CfD f UfVfDET 


A,B,Cf D,U, V 


- D* V ) /DET 

- C*U )/DET 


1 AUX0160 
IAUX0170 
1AUX01B0 
I AUX0190 
1AUX0200 
1AUX0210 
IAUX0220 
IAUX0230 
1 AUX0240 
1AUX0250 
1 AUX 0260 
IAUX0270 
1AUX0280 
IAUX0290 
1 AUX0300 
1AUX0310 
1 A JX0320 
i AUX0330 
1 AUX0340 
1AUX0350 
1AUX0360 
l AUX0370 
IAUX0380 
1AUX0390 
1 AUX0400 
1AUX0410 
l AUX0420 
1AUX0430 
1AUX0440 
1 AUX0450 
1AUX0460 
LAUX0470 
1AUX0480 
1 AUX0490 
1 A JX0500 
LDER0010 
1DER0020 
1DER0030 
1 DER0040 
1DER0050 
IDER0060 
1DER0070 
1DER0080 
i DER0090 
1DER0100 
1DER0110 
1DER01 20 
1DER0130 
1DER0140 
1DER01 50 
1DER0160 
1DER0 170 
1 DEROi 80 
1DER0190 
2AUX00I0 



n o o 


00 

to 


C MOD. 7/07/69 - L.F.H. 

C 

SUBROUTINE AUX2 

...DATA UPDATER FCR THE 1ST PHASE OUTPUT... 


C 


C 


c 


c 


c 


c 


c 


REALMS 



L (250) , 

1 



DUO) , M, JZ ( 3 ) , 

2 



A X 1 » Y ( 8 ) , DER ( 8 ) , 

3 



DC ( 5 ) » OS ( 5 ) * 

4 



F ( 6 ) » F2 3 * 

5 



FX( 3), FY(3) , 

6 



M 7 ( 3 ) , MY ( 3 ) » 

7 



CY ( 3 ) » CX ( 3 ) t 

8 



GAM r G AMD, TWPI r CRT 0 , C DTR , P l D2 , PID4 , 111 

COMMON 

/ 

DATA / 

Z»NEQL#NEQ2fJl,J?»J3tJ4,J5 

EQUIVALENCE 


I 



( 0(1), ZU1) ),( M , Z ( 2 1 ) ),( JZ ( 1 ) , Z ( 28 ) ), 

2 



( AX 1 1 Z ( 45 ) ),( Y ( 1 ) , Z ( 64 ) ),( DER(1),Z(9L) 

3 



( DC ( l ) , Z ( 13 l ) ), ( DS ( I ) , Z ( 1 36 ) ), 

3 



( F ( 1 ) , Z ( 141 ) ),( F 2 3 , Z ( 147) ), 

4 



( FX ( 1) , Z (151) ) , ( FY(1) , Z ( 154) ) f 

5 



( MZ(1),Z(161) ) t ( MY (1 ) , Z ( 1 64 ) ), 

6 



( CY(1),Z(171) ),( CXU),ZU74) ) , 

7 



( GAM , Z ( 32 ) ) ,( GAMD , Z ( 33 ) ) 

common 

/ 

CONS / 

TWPI,CRT0fCDTRfPID2,PID4 


F ( l ) = M*OER(2)*D(5) 

F(4) = M*DER(l)*D(7) 

F(5) = M*( Yll)**2 )*D(7) 


DC ( 3 ) - ( DU) - 0(5)*DC(1) )/D(7) 
DS ( 3 ) = D(5)*DSU)/n( 7) 


CALL F AX 1 
FXIl) = M# AX 1 

FX ( ? ) = FU)*DS<L) * F 2 3*DC ( L ) - F(4)*DS(3) 

L - F ( 5 ) *QC ( 3 ) - F(6)*DS<2) 

F X ( 3 ) = F(1)*DC(1) - F 2 3*DS ( 1 ) + F(4)*UC(3) 

L - F( 5) *D$( 3) - F (6) *DC(2) 

FY ( l ) = FX ( 1 ) 

F Y ( 2 ) = FX(2)*DC(2) - FX(3)*DS(2) 

F Y ( 3 ) = FX(3)*DC(2) + FX(2)*DS(2) 

111 = DFR ( I) - DER t 2) 

MZ( L ) = JZ( I)*ZZZ*DC (4 ) 

MZ(?) = ( J Z ( 3 ) - JZU) )*( ( Y(l) - Y { 2 ) ) **2 ) *DC ( A ) *DS l 4 ) 

MZ(3) = - JZ(3)*ZZZ*DS(4) 

MY ( 1 ) = MZ(1)*UC(4) - MZ(3)*DS(4) 

MY ( 2 ) = MZ(2) 

MY ( 3 ) = MZ(3)*0C(4) + MZU)*DS(4) 

CYtl) = MY ( l ) - F Y ( 2 ) *0 ( 2 ) - FY(3)*D<1) 


2AUX0020 
2AUX0030 
2AUX0040 
2 A JX 0 050 
2 AUX0060 
2 AUX0070 
2AUX0080 
2AUX0090 
2AUXOIOO 
2AUX0 l 10 
2AUX0120 
2AUX0130 
2AUX0140 
2AUX0150 
2AUX0160 
2AUX0170 
2AUX0180 
2AUX0190 
2 AUX0200 
2AUX0210 
2 AUX0220 
2 AUX0230 
2AUX0240 
2 AUXO 250 
2AUX0260 
2AUX0270 
2AUX0280 
2AUX0290 
2AUX0300 
2 AUX 0 310 
2 AUX0320 
2 AUX0330 
2 A'JXO 3 40 
2 AUX0350 
2 AUX0360 
2AJX0370 
2 AUX0380 
2 AUXO 390 
2 AUX0400 
2AUX0410 
2AUX0420 
2 AUX 0430 
2 AUX0440 
2AUX0450 
2 A JXQ 4 60 
2 AUX0470 
2 AUX0480 
2 AUX0490 
2 AUX0500 
2AUX0510 
2AJX052O 
2 AUX0530 
2 AUX0540 
2 AUX0550 
2AUX0560 



non n n n on oooo non on 


U> 

U) 


C Y( 2 ) = M Y( 2 ) + FY(1)*0(2) - FY(3)*0<3) + F(6)*D(4) 
CY ( 3 ) = MY t 3 ) v FY(1)*D<1) * FY(2>*D<3) 

C 

CX(1) = CYU) 

CX ( 2 ) -= CY ( 2 ) ♦DC ( 2 ) + CY(3)*DS(2) 

CX ( 3 ) = C Y ( 3 ) *0C ( 2 ) - C Y ( 2 ) *DS { 2 ) 

C 

GAMD = GAMX'CRTD 
C 

990 RETURN 
END 

2/18/69 - L.F.H. 

SUBROUTINE F AX L 

•..DUMMY LONG l T • ACCELERATION SU3P. . • 


RETURN 

END 

2/18/69 - L.F.H. 

MOD. 7/07/69 - L.F.H. 


SUBROUTINE FF6 


REAL*8 

X 

COMMON / DATA / 
EQUIVALENCE 
X 


2 ( 250 ) » 

F ( 6 J 

7^NEQl»NEQ2fJlfJ2*J3fJ4 t J5 
( F ( l ) , 7 ( 141 ) ) 


300 F ( 6 ) = 0.000 


990 RETURN 
END 

3/26/69 - L.F.H. 

MOD. 9/24/69 - L.F.H. 


SUBROUTINE OUT L ( JJ ) 

...OUTPUT SUBP. FCR THE 1ST PHASE... 

110 FORMAT ( 1 l* / »0», 39X, '...BASE DATA...* / IX 
X / 5X,‘T0 - INITIAL TIME ( SEC ).* 

X / 5 X * * HO - INITIAL INTEGRATOR TIME STEP SIZE ( SEC ).» 

X / 5X, * Q{ l ) . . . Q ( 3 ) - INITIAL OMEGA, ETA-00 T, AND ETA 1 

X * ( RPS, RPS, DEG ) .* 

X / 5X,'D( 1)...0(8) - CLAMSHELL GEOMETRIC PARAMETERS ( FT 
X / 5X,»M - CLAMSHELL MASS ( SLUG ).* 

X / 5X» 1 JC Y I D...JCY16) - CLAMSHELL Y-FRAME * 

X * MOMENTS AND PRODUCTS OF INERTIA ( SLUG*FT**2 ).» 

X / 5X,'JZ(l)...J/(3) - CLAMSHELL Z-FRAME I.E. PRINCIPAL * 
X ' MOMENTS OF INERTIA ( SLU0*FT**2 ).' 

X / 5X , * JVX t - CENTRAL BCDY ROLL MOMENT CF INERTIA * 

X •( $LUG*FT**2 ).' 

X / 5 X , * Y T - TERMINAL GAMMA ( DEG ).* 


2 AUX0570 
2AUX0580 
2AUX0590 
2 AUX0600 
2AUX0610 
2AUX0620 
2AUXD630 
2AUX0640 
2AUX0650 
2AUX0660 
2AUX0670 
FAXOOOIO 
FAX00020 
F AX00030 
FAX00040 
FAX00050 
FAX00060 
FAX00070 
FAX00080 
6FF00010 
6FF00020 
6FF00030 
6FF00040 
6FF00050 
6FF00060 
6FF00070 
6FF00080 
6FF00090 
6FF00100 
6FF00110 
6FF00120 
6FF001 30 
6FF00140 
6FF00150 
10UT0010 
10‘JT 0020 
IQUT0030 
TOUT 0040 
10UT0050 
10UT0060 
10UT0070 
1 0UT0080 
1 0UTQ090 
10UT01 00 
1QUT0U0 
1 GUT 0120 
10UT01 30 
IOUT0140 
10UT01 50 
I0UT0160 
iGUT 0170 
10UT01 80 
10UT0190 
IUUT0200 
IOUT0210 



o o 


4 ^ 


X / 5X f *YEPS - TERMINAL GAMMA CONVERGENCE C T I TER I ON ( DEG ).' 

X f 5X f 1 A X 1 - SYSTEM XI ACCELERATION ( FT/SEC**2 ).■ ) 

L 1 5 FORMAT! 1 1 1 / • 0* , 39X , • . . . BA SE DA TA • • . ' ) 

120 FORMAT! 'OCASE NO. ',I3,5X,10A8 / LX 
X / 6X» 1 TOfHOfO(J) • . • * »1P5012.4 

X / 12X,'Dl J)...* i 1 P 8D 1 2 • 4 

X / 8Xf*MtJCY(J) . • • ' tlP7D12.A 

X / 6X«'JZ(J) ,JVX1... # f lP4D12.4 

X / SX^YTiYEPStAXL... 1 ,1P3D12.4 ) 

130 FORMAT! • 1 1 / *OCA$E NO. ',I3 V 5X»L0A& 

X / 38X«'1ST PHASE OUTPUT' / IX 
X / 5Xt * T - ELAPSED TIME ( SEC 

X / 5X, 1 GAMD - CLAMSHELL ROLL-OUT ANGLE ( DEG ).' 

X / 5Xt * Y( 3 1 - CLAMSHELL ETA-ANGLE ( DEG ).' 

X / 5Xf 1 Y l 2 ) - CLAMSHELL ROLL-OUT RATE { RPS ).’ 

X / 5X* 1 Y ( L ) - SYSTEM ROLL RATE ( RPS ).' 

X / 5X, * DER ( 2 ) - CLAMSHELL ROLL-OUT ACCELERATION ( RAD/SEC++2 ). 

X / 5X f ' DER ( 1 ) - SYSTEM ROLL ACCELERATION ( RAD/SEC**2 ).* 

X / 5 X f 1 F ( J) - INERTIAL AND APPLIED FORCES ! LB 

X / 5Xf 1 F X ( J ) - X-FRAME HINGE FORCE COMPONENTS ( LB ),« 

X / 5X,'CX!J) - X-FRAME HINGE COUPLE COMPONENTS ( FT-L8 ).* 

X / 'O' / 55X , 1 . • . 1 

X / ' 0* , L5X, 1 T* ,9X,'GAMD' ,BX,* Y (3) 1 , 

X 8Xt , Yl2) i tBXf’YIl) 1 tYXt'DER^) 1 t6X f , DER(l) 1 

X / * 0* ,25X,' F( II • f 8X,' F ( 2 ) 1 ,8X, 1 F (3) • , 

X 8X » 1 F ( 4 ) 1 f 8X , f F { 5 ) * , 8X f • F (6 ) 1 

X / *0' 9 25Xf v FX(l)' «7X 9 v FX(2) v iTXt'FXOI* » 

X 7X« f CX( 1) ' ,7X,*CX<2) • ,7X,'CXI3) ' ) 

140 FORMAT! • L • / ■ OCASE NO. ■ ,I3,5X,10A8 

X / 38X i 1 1ST PHASE OUTPUT 1 ) 

150 FORMAT! IX / 10X,IP7012.4 / ( 22X,1P6D12.4 ) ) 


C 

C 


REAL+8 

X 

X 

X 

X 

X 

RE AL*8 
RE AL*8 

COMMON / CONS / 
COMMON / DATA / 
EQUIVALENCE 
X 
X 
X 
X 
X 
X 
X 
X 


2(250), 

TO f HO,Q! 3), D(8) ,M, JCYI6), 

JZ(3) • JVXl f GAMO, 

YT,YEPStAXl , 

XfYl3)fDER(3)|F(6),FX(3) f CX(3)t 
TITLE! 10) 

YY ( 3 ) 

TWPI,CRTDfCDTRfPl02fPID4 
TWPI t CRT0 f C0TR,PID2,PID4 
Z,NEOl ,NEQ2 , Jl f J2 , J3, J4 f J5 

( T 0, Z ! I ) ),( H0 V Z ( 2 ) ) f (.Q(i),Z(3) )• 

( D ( 1 ) « Z ( 1 1 ) ),( M v Z ( 2 1 ) ),( JCY ( 1 ) , Z ( 22 ) ), 

( J Z ( 1 ) , Z ( 28 ) ),( J VX 1 , Z ( 3 1 ) ) t ! G AMD , Z ! 33 ) ), 
( YT,Z!42) >,( YEPS,Z(43> ),( AXl f Z(45) ), 

( X , Z ( 6 1 ) ),( Y ( 1 ) , Z ( 64 ) ),( DER(l),Z(9l) ), 
l F ( 1 ) , Z ( 141 ) ), 

( FX< 1) ,Z( 151) ) , ( CX( 1), Z( 174) ) f 
( T ITLE ( 1 ) , Z (241 ) ) 


DATA 


JK/ 0 /, NCASE / 0 / 


10UT0220 

l 0UT0230 

1 0UT0240 

10UT0250 

1 0UTC260 

10UT 0270 

10UT0280 

10UT0290 

10UT0300 

1 OUT031 0 

10UT0320 

10UT0330 

10UT0340 

I0UT0350 

10UT0360 

10UT0370 

1 OUT 0380 

IOUT0390 

10UT 0400 

IOUT0410 

10UT0420 

10UT0430 

10UT0440 

10UT 0450 

10UT0460 

10UT 0470 

1 OUT 0480 

10UT0490 

10UT0500 

1QUT051 0 

10UT 0520 

l 0UT0530 

1 OUT0540 

10UT0550 

10UT0560 

10UT0570 

IOUT0580 

10UT05 90 

1 0UT0600 

1 0UT061 0 

10UT0620 

10UT0630 

10UT 0640 

10UT0650 

1 0UT0660 

10UT0670 

1DUT0680 

10UT0690 

10UT0700 

10UT0710 

10UT0720 

10UT0730 

10UT0740 

1 0UT0750 

1QUT0760 



non 


I F { JJ . NE. 0 ) 

GO TJ 400 

10UT0770 

IFI NC AS E .NE. 

J 4 ) GO TO 230 

IOJT0780 

NCASE = J4 + i 


10UT0790 

J4 = NCASE 


1 OUT0800 

GO TO 240 


10UT0810 

230 NCASE = J4 


1 0UT0820 

240 K F ( JK .EQ. 0 ) 

WRITE! 6,110 ) 

10UT0830 

250 I F ( JK .NE. 0 1 

WRITE! 6,115 ) 

10UT0840 

10UT0850 

WRITE! 6,120 ) 

NCASE, TITLE, 

10UT 0860 

X 

T 0 , HO , Q, 

10UT0870 

X 

D, 

1 0UTO88O 

X 

M , JCY , 

10UT0890 

X 

J7, JVXlt 

1 0UT0900 

X 

YT , YCPS, AX 1 

1 OUT0910 

300 JK = 255 


10UT092 0 
10UT0930 

WRITE! 6,130 ) 

NCASE, TITLE 

10UT0940 

JL = 20 


1 OUT 0950 
10UT0960 

400 I F ( MOO ( JL , 56 

) .NE. 0 ) GO TO 500 

1 QUT0970 

WRITE! 6,140 ) 

NCASE , TITLE 

1 0UT0980 

JL = 0 


10UT0990 
10UT 1000 

500 YYI 1 ) = YID/TWPI 

10UT1010 

YY ( 2 J = Y ( 2 ) /TWP I 

1QUT1020 

YY ( 3 ) = Y ( 3 ) *CRTD 

l OUT 1030 

WR 1 TE { 6,150 ) 


10UT1040 

X 

X,GAMD,YY(3) ,YY!2) ,YY( 1),DER(2> ,DER(1 ) , 

10 JT 1 050 

X 

F, 

10UT1060 

X 

FX,CX 

10UT 1 070 

JL = JL + 4 


10UT 1080 
10UT1090 

990 RETURN 


10UTL100 

END 


10UT1110 

2/26/69 - L.F.H. 


2NIT0010 

MOD. 9/03/69 - L.F. 

H. 

?N I T 0020 
2NIT0030 

SUBROUTINE NIT2 


2NIT0040 

2NIT0050 

...DATA 

INIT. FOR THE 2ND PHASE... 

2NIT0060 

2NIT0Q70 

REAL*8 

/ ( 250 ) , 

2N I TOO 80 

X 

DllO) , 

2N I T 0090 

X 

GAM, XT, 

2NIT0100 

X 

H, BETA, 

2NIT0110 

X 

XF, 

2NIT0120 

X 

X,Y{8),W(8>, 

2N I TO 1 30 

X 

DC ( 5 ) , DS ( 5 ) , 

2N1T0L4O 

X 

D7CA , D7S A , 

2NIT0150 

X 

DXCM ( 3 ) , P ( 3, 4 ) , 

2NIT0160 

X 

DaM 3 ,DtM 2 ,YiM 2 ,ZZZ, 

2NIT0170 

X 

TWPI ,CRTD,CDTR,PID2,PID4 

2NIT0180 

COMMON / DATA / 

Z,NEQ1,NEQ2,J1,J2,J3,J4,J5 

2NIT0190 

equivalence 


2NIT0200 



u> 

as 


x I D( 1),ZI m ) , 

X ( G AM f Z ( 32 ) ),( XT * Z ( 41 ) ), 

X ( H, 7.144) ),( BETA,Z(47) ), 

X ( XF.ZISl) ), 

X ( X, 7(61) ), ( Y( 1 ) # W( 1) ,7(64) ) , 

X ( DCU), 7.(131) ),( QS ( 1 ) » Z ( 1 36 ) )* 

X ( 0 7C A * Z ( 14 1 ) ),( 07SA f Z(142) ), 

X ( 0XCH(l)|Z(t51) ) f ( P(l»l)*Z(16l ) ) 

C 

COMMON / CONS / T WP t , CRT D f CD I R , P I D 2 , P I D4 
C 

EQUIVALENCE l 03M 3 f D IM2 , Y1M2 , l Z Z ) 

C 

XF = X 
C 

C « • • BODY POINT DISPLACEMENTS IN THE Z-FRAME... 

C 

D8M3 = DCS) - D ( 3 ) 

P(ltl) = D8M3* DC ( 4 ) - D(2)*DS(4) 

P ( 2 » 1 ) = O.ODO 

P ( 3 » 1 ) = - D(2)*DC(4) - 08M 3*0$ ( 4 ) 

D 1 M2 = D(l) - D( 2) 

P ( 1 1 2 ) = - D(3)*DC(4) +■ D 1 K2*D S ( 4 ) 

P ( 2 f 2 ) = O.ODO 

P ( 3 * 2 ) = DI M2*GC ( 4 ) + D ( 3 ) *DS ( 4 ) 

PI 1*3) = - D(3)*DC(4> - D(2)*DS(4> 

P ( 2 » 3 ) = Dll) 

P ( 3 i 3) = - D ( 2 ) * DC ( 4 ) + D ( 3 ) *0 S ( 4 ) 

P( L » 4 ) = P( l t 3) 

P 1 2 f 4 ) = - Dll) 

P ( 3*4) = P ( 3*3) 

C 

C ...BODY C.M. DISPLACEMENT RATES... 

C 

D7C A = D(7)*DC(3) 

D7SA = D( 7) *DS( 3) 

DXC M ( 2 ) = Y(2)*D(5)*DS(l) - Y(L)*D7SA 
DXCMI3) = Y(2)*D(5)*0C<1) + Y(t)*D7CA 
C 

C ...INIT. BODY FIXED Z-FRAMF ANGULAR RATES... 

C 

YIM2 = Y( t) - Y ( 2 ) 

W(l) = Y1M2*DC(4) 

W t 2 ) = O.ODO 
W ( 3 ) = - Y1M2*DS(4) 

C 

C ...INIT. EULER I AN ANGLES... 

C SEE NOTES IN SUBP. UER2 AND SUBP. AUX4 ON ROTATIONAL TRAN SF OR MA T I ON. 
C 

W ( 4 ) = - GAM 
W(5> = O.ODO 
W( 6) = - BETA 
C 

C ...SET UP TIME STOP AND ' H 1 . . . 

C 


2NIT0210 

2NIT0220 

2NIT0230 

2NIT0240 

2NIT0250 

2NIT0260 

2NIT0270 

2NITO20O 

2NIT0290 

2NI T0300 

2NIT0310 

2N I TO 320 

2NIT0330 

2N I T 0340 

2NIT0350 

2 MI TO 360 

2N l TO 370 

2NIT0380 

2NIT0390 

2NIT0400 

2NIT0410 

2NIT0420 

2N I T 0430 

2NIT0440 

2NIT0450 

2NIT0460 

2NIT0470 

2NITQ480 

2NIT0490 

2NIT0500 

2NIT0510 

2NIT0520 

2N I T 0530 

2NIT0540 

2NITU550 

2NIT0560 

2NIT0570 

2NIT0580 

2NIT0590 

2NIT0600 

2NIT0610 

2NI T0620 

2NIT0630 

2NIT0640 

2NIT0650 

2NIT0660 

2NIT0670 

2NIT0680 

2NIT0690 

2N I TO 700 

2N I TO 710 

2NIT0720 

2NIT0730 

2NIT0740 

2NIT0750 



<Js) 

<1 

I 


III - DA BS ( Y 1 M2 ) 

I F { III . LT , PIU4 ) III =■ PI04 
III = PI 02/ III 
XT = X + 2£Z 
H * ZZZ/128,000 
C 

990 RETURN 
FND 

C 2/2T/69 - L.F.H. 

C MOD, 8/20/69 - L.F.H. 

C 

SUBROUTINE $ET2( T , W , DER f NCO , H t J 1 ) 

C 

C ...TERMINAL H CONTROLLER FOR THE 2 4 D PHASE... 

C 

REAL*8 ZI250), 

X TT»H f 

X TTNeO) ,W(NEQ,3) ,OFR(NEO,A) , 

X Hl,H2 

COMMON / DATA / Z,KKK(6) f J5 
EQUIVALENCE 

X ( T T » Z ( 4 1 ) ) 

C 

IF! J1 .EQ. 0 ) JC = 0 
C 

I F { J1 .EQ. 3 ) GO TO 255 
IF ( J5 .EQ. 0 ) GO TO 250 
J5 = 0 

CALL SPIEl 0 ) 

JC = JC + I 

I F ( JC .GE. 5 ) GO TO 500 
HI = H 
H = 2.0D0*H 
T(l) ? T ( 3 ) 

DO 200 K = 1 » NFQ 
200 W ( Kt 1 ) = W ( K f 3 ) 

Jl = 3 
C 

230 I F { T ( I ) .EQ. TT ) GO TO 300 

IF { Til) + H .LE. TT ) GO TO 990 
H2 = TT - T< 1) 

JC = 0 

I F ( Jl .EQ. 3 .AND. HI .EG. H2 ) GO TO SOO 
H = H2 
GO TO 990 
C 

250 I F ( JC .EG. 0 ) GO TO 255 
H = H*( 0.5D0**JC ) 

JC = 0 
255 Jl = 1 

GO TO 230 
C 

300 Jl = 2 

GO TO 990 
C 


2NIT076Q 
2NIT0770 
2NIT0780 
2N I T 0790 
2NIT0800 
2NIT0810 
2NIT0820 
2NIT0830 
2SET00 10 
2SET0020 
2SET0030 
2SET 0040 
2SET0050 
2SFT0060 
2SET0070 
2SET00 80 
2SET0090 
2SET0100 
2SET0110 
2SET0120 
2SET0130 
2SET0L4O 
2SET0150 
2SET0160 
2 SE TO l 70 
2SET0180 
2 SE TO 190 
2SET0200 
2SET021 0 
2SFT0220 
2SET0230 
2SET0240 
2SET0250 
2SFTD260 
2SET0270 
2SET0280 
2SETU290 
2SET0300 
2 SFT0310 
2SET0320 
2SET0330 
2SET034G 
2SET0350 
2SET0360 
2SET0370 
2SET0380 
2 SE r 0390 
2SET0400 
2SET0410 
2SET0420 
2SET0430 
2SET0440 
2SET0450 
2SE TO 460 
2SET0470 



u> 

00 


500 Jl = 5 
990 RETURN 
END 

C 2/26/69 - L.F.H. 

C MOD. 7/07/69 - L.F.H. 
C 

SUBROUTINE AUX3 


C 

c 

c 


c 

c 

c 


c 

c 

c 


c 


...DATA UPDATER FOR THE 2ND PHASE DERIV. SUBP... 


REAL*8 


COMMON / DATA / 
EQUIVALENCE 


Z. ( 2 5 0 ) f 

W ( 8 ) » 

DC I 5 ) » DS ( 5) 

Z,NE01,NEQ2 f Jl* J2, J3 f J4.J5 
{ Wtn f Z(66) ), 

I DC ( l ) t Z ( 13 L ) ) w ( DS (1 ) * Z 1 1 36 ) ) 


...COSINES AND SINES OF THE EULFRIAN ANGLES... 


DO 300 J = 1,3 

DC I J ) = DCOS ( ( J + 3 ) ) 

300 DS(J) = DSINl WtJ+3) ) 

...CHECK THE COSINE OF THE 2-ANGLE... 


1 F ( DC ( 2 ) .NF. O.ODO ) GO TO 990 
CALL SPIEI 1 ) 

J 5 = 255 


990 RETURN 
END 

C 2/24/69 - L.F.H. 

C MOD. 7/07/69 - L.F.H. 
C 


C 

C 

C 


C 

C 

c 


SUBROUTINE DER2 { T » W , DER ,NEQ , J ) 


...DERIV. EQ'S. FOR THE 2ND PHASE... 


RE AL *8 
REAL *8 

1 

COMMON / DATA / 
EQUIVALENCE 

1 

2 


T » W ( NEQ ) fDER(NEQ,4) 

Z ( 250 ) t 

J7(3)#DC(5)fDS(5) 

Z,NEQ1,NEQ2, Jl, J2 f J3, JA t J5 

( JZ I 1 ) , Z (28) ),( DCfl)«ZU31) )• 
( DS(l).Z<136) ) 


...EULER'S EQ'S. OF MOTION FOR FREE FLIGHT... 


DERI 1 , J ) = W t 2 ) *W ( 3 ) * I J Z ( 2 ) - JZ(3) )/JZ(I) 

DER ( 2, J ) = W(3)*W<1)*( J Z ( 3 ) - JZ I 1 ) )/JZ(2) 

DER I 3 » J ) = Hfl)*W(2)*( JZ(1) - J Z ( 2 ) J/JZI3) 

C 

C ...DERIV. EQ'S. FOR EULERI AN ANGLES... 

C THE ROTATIONAL TRANSFORMATION FROM THE INERTIAL FRAME TO THE 


2SET0480 

2 SET 0 A 90 

2SET0500 

3AUX0010 

3AUX0020 

3AUX0030 

3 AUX0040 

3AUX0050 

3AUX00S0 

3AUX0070 

3AUX0080 

3 AUX0090 

3AUX0100 

3AUX0110 

3AUX0120 

3AUXOL30 

3AUX0140 

3AUX0150 

3AUX0160 

3AUX0170 

3AUX0 180 

3AUX0190 

3AUX0200 

3AUX021 0 

3 AJX0220 

3AUX0230 

3AUX0240 

3 AUX 0250 

3AUX0260 

3AUX0270 

3AUX0280 

3AUX0290 

2DER0010 

2DER0020 

2DER0030 

2DER0040 

2DER0050 

2DER0060 

2DER0070 

2DER0080 

2DER0090 

2DER0100 

2DER0110 

2DER0120 

2DER0130 

2DER0140 

2DER0150 

2DER0160 

2DER01 70 

2DER01 80 

2DER0190 

2 DER0200 

2DER0210 

2DER0220 

2DER0230 



nnoooooooo ooo o r> o 


U) 

VO 


BODY FIXED PRINCIPAL 
SUCCESSIVE ROTATIONS: 
A ROTATION ABOUT THT 
A ROTATION ABOUT THF 


AXIS FRAME AT THF BODY C.M. INVOLVES THREE 
A ROTATION ABOUT THF 1-AXIS THROUGH THE I- 
1ST INTERMEDIATE 3-AXIS THROUGH THE 2-ANGLE 
2ND INTERMEDIATE 2-AXIS THROUGH THE 3-ANGLE 


) /DC ( 2 ) 


DER ( 4, J ) = ( W I 3 ) *D S ( 3 ) + W(1)*DC13) 
DER ( 5, J ) = W(3)*DCU) - W(l)*DS<3) 
DER ( 6 f J ) = W 1 2 T + DER ( 4, J ) *D$( 2 ) 

990 RETURN 
END 

2/24/69 - L.F.H. 

MOD. 9/03/69 - L.F.H. 


SUBROUTINE AUX4 

...DATA UPDATER FOR THE 2ND PHASE UUTPJT... 

REAL+8 Z (250) v 

X TF.TD, 

X T f W ( 0 ) f 

X DC ( 5 ) , D$( 5) , 

X D7CA,07SA, 

X 0 XC M ( 3 ) , XC M ( 3 > , 

X P(3,4) ,X(3,4) » 

X A I 3f 3 ) 

Z , NEC 1 , N 1 02 , J 1 , J? « J 3 , J4, J 5 


COMMON / DATA 
EQUIVALENCE 


( TF f Z ( 5 l ) ), 

( T , Z ( 6 1 ) ), ( W ( 1 ) v Z(64) ), 

t DCIl)»ZU3l) ),l DS ( 1 ) , 7 I 1 36 ) ), 

( D7CA,Z( 1411 I f ( D 7S A , Z ( 142) ), 

( DXCM(1)«Z(L51) >,( XCM ( 1 ) , Z ( l 54 ) ), 
( P(l«l)«Z(16l) >,< X ( 1 « I ) * Z ( ?01 ) ) 


...BODY C.M. DISPLACEMENTS... 

XCM ( 1 ) = O.ODO 
TD = T - TF 

XCM ( 2 ) = DXCM(2)*TD + D7C A 
XCM t 3) = OXCM(3)*TO + D7SA 

THE ROTATIONAL TRANSFORMATION FROM THE INERTIAL FRAME TO THE 
BODY FIXED PRINCIPAL AXIS FRAME AT THE BUDY C.M. INVOLVES THREE 
SUCCESSIVE ROTATIONS: A ROTATION ABOUT THE 1-AXIS THROUGH THE 1- 
A ROTATION ABOUT THE 1 ST INTERMEDIATE 3-AXIS THROUGH THE 2-ANGLE 
A ROTATION ABOUT THE 2ND INTERMEDIATE 2-AXIS THROUGH THE 3-ANSLE 
THE FOLLOWING ELEMENTS ARE COMPONENTS OF THF INVERSE OF THE 
FORMATION MATRIX FOR ROTATING THE INERTIAL FRAME INTO THE BODY F 
PRINCIPAL AXIS FRAME AT THE C.M. OF THE BUDY. 

A ( 1 9 1 ) = DC(2)*DCm 
A ( l , 2 ) = - DS ( 2 ) 

A ( 1,3) = DC ( 2 ) *DS I 3 ) 

A ( 2 , 1 ) = DStl)* 0 St 3 ) + DCm*DS(?)*QC(3) 


2DER0240 
ANGLE, 2DER0250 
AND 2DER0260 
2DER0270 
2DER0280 
2DFR0290 
2DER0300 
2DER0310 
2DER0320 
2DER0330 
2DER0340 
4AUX0010 
4AUX0320 
4AUX0030 
4AUX0040 
4AUX0050 
4AUX0060 
4AUX0070 
4 AUX0080 
4AUX0090 
4AUX0100 
4AUX0110 
4AJX0120 
4AUX0130 
4AUX0140 
4AUX0IS0 
4AUX0160 
4AUX0170 
4AUX0 180 
4AUX01 90 
4AUX0200 
4AUX0210 
4AUX0220 
4AUX0230 
4 AUX0240 
4AUX0250 
4AUX0260 
4AUX02 70 
4AUX0280 
4 AUX0290 
4 AUX0300 
4AUX031 0 
4AUX0320 
4AUX0330 
ANGLE, 4AUX0340 
, AND 4AUX0350 
. 4AUX0360 

TRANS-4AUX0370 
IXhD 4 AUX0380 
4AUX0390 
4 AUX0400 
4AUX0410 
4AUX0420 
4AUX0430 
4 AUX 0440 



non non o ooo 


-P* 

O 


A (2, 2) = DC ( 1 ) *0C ( ? ) 

A ( 2 , 3 ) = - D S ( 1 ) +DC ( 3 ) + DC ( l ) *DS 12 ) *DS ( 3 ) 
A ( 3 , l ) = - DC ( 1 ) *DS f 3 ) + DS(1)*DS(2)*DC(3) 
A ( 3 , 2 ) = DS(1J*DC12) 

A ( 3,3) = DC ( 1 ) *DC ( 3 ) + US ( t ) *DS ( 2 ) *05 ( 3 ) 
•..BODY POINT DISPLACEMENTS... 


DO 500 J = 1,4 
DO 500 K = I. 3 
X ( K f J ) = XCM(K) 

DO 500 L = 1,3 

500 X ( K , J ) = X ( K , J ) + A(K,L)*P(L, J) 

990 RETURN 
END 

3/26/69 - l.F.H. 

MOD. 9/03/69 - L.F.H. 


SUBROUTINE 0UT2( JJ ) 


...OUTPUT SUBP. FOR THE 2ND PHASE... 


C 


C 


c 


L30 FORMAT ( 'l 1 / 'OCASE NO. ',I3,5X,10A8 

X / 38X , 1 2ND PHASE OUTPUT' / IX 

X / 5 X , 1 T - ELAPSED TIME ( SEC ).' 

X / 5X, ' XCM , YCM,2CM - C.M. FREE FLIGHT DISPLACEMENT ( FT ).' 

X / 5X,'Xl J) ,Y( J) ,2(J) - J-TH POINT FREE FLIGHT DISPLACEMENTS 

X M FT ).' 

X / 'O' / 55X , * . . . ' 

X / ' O' ,14X,' T' tlOXr'XCM' #9X , ' YCM' ,9X,'2CM', 

X 9X,'X(1)* ,8X,'Y(I)' ,8X,'Z(l)' 

X / 'O' ,25X,'X(2)' ,0X,'Y(2)' ,8X,'2<2)' , 

X SXf'XIB)' ,8X,'Y(3) # , 8 X , ' Z ( 3 I ' 

X / 'O' ,25X,'X(4)' ,8X,'Y(4)' ,8X,'Z(4)' ) 

140 FORMAT ( '1' / 'OCASE NO. ' , I 3 , 5X , 10 A8 

X / 3 8X , ' 2ND PHASE OUTPUT' ) 

150 FORMAT ( ' 0' , 9X , 1P7D 12 .4 / ( 22X , 1 P601 ? . 4 ) ) 


REAL*8 

X 

X 

X 

COMMON / DATA 
EQUIVALENCE 
X 
X 
X 
X 


2(250), 

T, 

XCM ( 3 ) , X ( 3, 4 ) , 

TITLE! 10) 

Z,NEGI,NE02,J1,J2,J3,J4,J5 
( T , Z ( 61 ) ), 

( XCM (!),/( 154) ),( XII, 1), 2(201) ), 
( TITLE* l), 2(241) ), 

( NCAS E , J4 ) 


I F ( JJ .NE. 0 ) GO TO 400 
WR I TE ( 6,130 ) NCASE, TITLE 
JL = 16 


400 I F ( MOD ( JL , 56 ) .NF. 0 > GO TO 500 


4 AUXO 450 

4AUX0460 

4 AUX0470 

4 AUX0480 

4AUX0490 

4 AUX0500 

4AUX05L0 

4AUX0520 

4AUX0530 

4AUX0540 

4AUX0550 

4AUX0560 

4AUX0570 

4AUX0580 

4AUX0590 

4AUX0600 

20UT0010 

20UT0020 

2QUT0030 

2QUT0040 

20UT0050 

20UT0060 

20 JT 0070 

20UT0080 

20UT0090 

20UT0100 

20UT0110 

20UT0120 

20UT01 30 

20UT0140 

20UT0150 

20UT01 60 

2UUT0L70 

2GUT0180 

20UT0190 

20UT0200 

20UT021 0 

2GUT0220 

20UT0230 

20UT0240 

20UT0250 

20UT0260 

20UT0270 

20UT0280 

20UT0290 

20UT0300 

20UT0310 

2 0UT0320 

20UT0330 

20UT0340 

20UT0350 

20UT0360 

20UT0370 

20UT0380 

20UT0390 



* * 


WR I TE ( 6,140 ) NCASE, TITLE 

JL = 0 
C 

500 WRITE! 6,150 ) T,XCH,X 

JL - JL + 4 
C 

990 RETURN 
END 

/* 

ft EXEC ASSEMBLR, PAR M=* LOAD, DECK* 

//SOURCE. SYSPUNCH DD DSN-CDEC K ,SYSOUT=B 
//SOURCE. SYSIN UD * 

EDAT TITLE 'EDAT A-ENT R Y AND INITIALIZATION SECTS.' 

PRINT DATA 

* 2/22/67 - L.F.H. 

* MOD. S/15/70 - L.F.H. 

SPACE 1 

* SUDP. EDATA PROCESSES DATA CARDS 'READ' BY SUBP. LOAD AND STORES 

* THE PROCESSED DATA ITEMS AS OIRECTCD. CHARACTER, FULL AND DOUBLE 

* WORD DECIMAL ( WHICH INCL. FULL WORD INTEGER AS WELL AS SHORT AND 

* LONG RFAL ), HALF WORD DECIMAL INTEGER, AND FULL WORD HEXADECIMAL 

* DATA ARE PROCESSED. 

* SUBP. EDATA IS MOST EFFICIENT WHEN DIRECTFO TOWARD THE LOADING OF 

* BLOCKS OF STORAGE SUCH AS A COMMON AREA OR A LARGE ARRAY CONTAINING 

* THE VARIOUS TYPES UF DATA REQUIRED BY A PRJCFSSING PROGRAM SINCE IT 

* TENDS TO LOAD THE PROCESSED DATA ITEMS INTO SEQUENTIALLY HIGHER AND 

* HIGHFR LOCATIONS. IT MAY BE DIRECTED TO SCATTER LOAD ONLY THE RE- 

* QUIREO DATA ITEMS, I.E., ALL OF THE STORAGE BLJCK NEED NOT BE LOADED. 

* THE CURRENT LOAD POINT MAY BE SHIFTFD AS RFG'D TO SKIP OVER SECTIONS 

* OF STORAGE AT ANY TIME. THIS SHIFT IS COMPUTED RELATIVE TO THE 
+ INITIAL LOAD POINT. 

* CHARACTER DATA ARE STORED WITHOUT CHANGE BYTF-3Y-BYTE ; HENCE, THE 

* CURRENT LOAD POINT CAN BE SHIFTED OFF A HALF WORD AS WELL AS A FULL 

* WURD AND DOUBLE WORD BOUNDARY. IT IS ADVISABLE TO USE THE INITIAL 

* LOAD POINT RESET AND SHIFT FEATURE OF SUBP. FDATA FOLLOWING THE LOAD- 

ING OF CHARACTER DATA. 

THE LOAD POINT FUR THE FULL AND DOUBLE WORD DECIMAL DATA AS WELL 

* AS THE HEXADECIMAL AND HALF WORD DECIMAL INTEGER DATA MUST BE AT THE 

* APPROPRIATE BYTE BOUNDARIES. 

* SUBP. FDATA PROCESSES ITS OWN PROGRAM INTERRUPTS. IT SAVES THE 

* ADDRESS OF 'THE NEXT INSTRUCTION* AND THE OLD PSW AND RESETS THE NEXT 

* INSTRUCTION ADDRESS BEFORE RETURNING CONTROL TO THE CONTROL PROGRAM. 

* THE CARD ERROR COLUMN NUMBER, THE CARD IMAGE, AND THE OLD PSW ARE 

* PRINTED BY SUBP. LOAD. 

* SUBP. EDATA CAN DETECT CERTAIN PROCESSING ERRORS AND GENERATE THE 

* CARD COLUMN NUMBFR AT WHICH THE ERROR WAS DETECTFD. THE LOADING OF 

* ALL SUBSEQUENT DATA IS SUSPENDED, AS IT IS WHEN A PROG. INTERRUPT IS 

* PROCESSED. HOWEVER, SUBSEQUENT DATA CAROS ARE CHECKED FOR ERRORS 

* UNTIL A RETURN TO THE CALLER OF SUBP. LOAD IS EXECUTED. ONLY ONE 

* ERROR PER CARD CAN BE DETECTED AND PROCESSED. DATA ITEMS ON A GIVEN 

* CARD WHICH FOLLOW A DETECTED ERROR CANNUT BE CHECKED. SUBP. LOAD 

* PRINTS THE CARD ERROR COLUMN NUMBFR AND THE CARD IMAGE IN SUCH CASES. 

* FOR DETAILS ON SUBP. LOAD AND SUBP. EDATA, CALL L .F .HAT AKEY AMA, 

* NASA-GSFC, CODE 721, X4047. 

SPACE 1 


20UT0400 
20UT0410 
20UT0420 
20UT0430 
20UT0440 
2 0UT0450 
20UT0460 
20UT0470 


EDATOOIO 
EDAT 0020 
EDAT0030 
ED AT 0040 
EDAT0050 
EDAT0060 
EDAT 0070 
EDAT0080 
EDAT 0090 
EDATOLOO 
EDAT0110 
EDAT0120 
EDAT0130 
EDAT0140 
ED A 10150 
EDAT0160 
EDAT0170 
EDAT01 80 
E0AT0190 
EDA T0200 
EDAT0210 
EDAT0220 
EDAT0230 
EDAT0240 
EDAT0250 
EDAT0260 
EDAT0270 
EDAT0280 
EDAT0290 
EDAT 0300 
FDAT0310 
EDAT 0320 
EDATO 330 
EDAT0340 
EDAT 0350 
EDAT 0360 
EDAT0370 
EDAT0380 
EDAT0390 
EDAT0400 
EDAT041 0 
EDAT0420 
EDAT 0430 



EDATA 

CSECT 



E0AT0440 


SPACE 

I 


EDAT0450 

A 

EQU 

10 


EDAT0460 

B 

FQU 

ll 


EDAT0470 

C 

EQU 

12 


EDAT0480 

0 

EQU 

13 


EDAT0490 

E 

FQU 

14 


EDAT0500 

F 

EQU 

15 


EDAT051 0 


EJECT 



EDAT0520 

* 


* * * 


ED AT05.30 


ENTRY 

BUST 


FDAT0540 


USING 

F 


FDAT0550 


R 

A100 


E DAT 0560 


OC 

CL6* 5EDATA* 


EDAT 0570 

ALOO 

STM 

E#C»12(D) 

SAVE THE R * S . 

EDAT0580 


L 

B, 4 ( D) 

GET THE PREV. SAVE AREA ADDR... 

EDAT 0590 


L 

A, 24(B) 

...C THE PARM. LIST ADDR. 

EDAT0600 


L 

4 , =A ( LFH2 ) 

SET UP THE CURRENT SAVE AREA... 

EDAT 061 0 


ST 

D * 4 ( 4 ) 

• • • 

EDAT0620 


ST 

4 1 8 ( D ) 

• • • 

EDAT0630 


LR 

0*4 

• • • 

EDAT 0640 


SPACE 

1 


EDAT0650 


LA 

F v BUST-EDAT A ( F ) 

SHIFT THE BASE ADDRESS... 

EDAT 0 660 


USING 

BUST « F 

• • • 

E0AT0670 


B 

A210 


EDAT0680 

♦ 


* * 

******** 

EDAT 0690 

OUST 

B 

A200 


EDAT0700 


DC 

CL6 1 5BUST * 


EDAT0710 

A200 

N I 

A212+1 ,X* OF* 

DISABLE THE BR. AT A212. 

EDAT0720 


SPACE 

l 


EDAT0730 

A210 

LR 

2,F 

SHIFT BASE REGISTERS... 

EDAT0740 


USING 

BUST, 2 

« * • 

EDAT0750 


DRUP 

F 


EDAT0760 


SPACE 

1 


EDAT0770 

A212 

RC 

15, A300 

** 

EDAT0780 


SPACE 

1 


EDAT0790 


MVC 

176(8, 4), 4(1) 

SET UP THE O-PSW. 

FDAT0800 


MVC 

9(3,1), =AL3(C110) 

MOD. THE O-PSW NEXT INSTR. ADDR. 

EDATOBIO 


01 

A212+l,X*F0* 

ENABLE THE BR. AT A212. 

EDAT0820 


BCR 

15, E 

RFTURN TO THE CONTOL PROG. 

EDAT0830 

* 


* * 

******** 

EDAT0840 

A300 

SPIE 

BUST, (11,15)) 

SPEC.. A NEW P. I.E. 

EDAT0850 


SPACE 

1 


EDAT0860 


STM 

1,2,24(4) 

SAVE R 1 C R2. 

EDAT0870 


CLC 

82(2,4>,=H*1* 

CHECK THE CARD COUNT. 

EDAT0880 


BH 

A400 

...BR. OUT IF NOT THE 1ST CARD. 

EDAT0890 


L 

3, 0(A) 

GET THE INI T . LOAD POINT ADDR... 

EDAT0900 


ST 

3,172(4) 

...STASH IT FOR FUTURE REF. 

EDAT0910 


LA 

5,88(4) 

LOAD R5 WITH A(CARD). 

EDAT0920 


LA 

7,79(5) 

SET UP R 7 WITH A ( CARD + 79 ) . • • 

EOAT0930 


LA 

8, LOB 

...RB WITH A ( LOB ) 

EDAT 0940 


LA 

9,3 

.. .C R9 WITH 3. 

EDAT0950 


STM 

1,9,24(4) 

SAVE R1 - R9. 

EDAT0960 


EJECT 



EDAT0970 

A400 

LM 

1,9,24(4) 

RESET R1 - R9. 

EDAT0930 




LR 

6,5 

RESET £ SAVE R6... 

EDAT0990 


ST 

6 ,44(4) 

• • # 

EDAT1000 


XC 

17618, A), 176(A) 

CLEAR THE *PSW* CELLS. 

EDAT 1010 


xc 

CL 8 ,C L 8 

ZERO LOP £ HOLD. 

EDAT1020 


LA 

A, 3 

SET UP R A t RB, RE, £ RF... 

EDATl 030 


LA 

B, A402 

... 

EDAT10A0 


LM 

E , F , LOP 

... 

FDAT1050 


01 

A408+ l , X * FO * 

ENABLE THE OR. AT AA08. 

EDATl 060 


SPACF 

L 


EDATl 070 

A402 

CLI 

0(6), C* * 

COMPUTE THF LOADING OFFSET... 

EDA T 1 080 


RE 

AA07 

• • • 

fcDATl 090 


NI 

AA08+ I , X 1 OF 1 

DISABLE THE BR. AT AA08. 

EDATl LOO 


LA 

B.AAOA 

...RESET RB. 

EDATl HO 

A404 

CLI 

0(6), C'O' 

« # « 

EDATl 120 


BL 

AA 05 

• • • 

EDAT 1130 


CLI 

0(6) ,C*9* 

• • • 

EDATl 140 


RH 

CIIO 

• • « 

EDATl 1 50 


MVN 

0(1, 8), 0(6) 

• mm 

EDAT 1 l 60 


M 

E , =F * 10* 

m m m 

EDATl 17 0 


A 

F , HOLD 

m • • 

EDATl 180 


B 

AA07 


F DA T 1 1 90 

A405 

LA 

B , AA 06 

...RESET RB. 

EDAT 1200 

A406 

CLI 

016) ,C • » 

• mm 

EDATl 21 0 


BNE 

CILO 

• • • 

EDATl 220 

A407 

LA 

6,1(6) 

• mm 

EDAT 1230 


RCTR 

A , B 


ED AT 1 240 


SPACE 

1 


EDATl 250 


CLI 

0(61,0' 1 

CHECK THF COL. 9 CHAR. 

EDAT l 260 


BNE 

Clio 


EDATl 270 


SPACE 

1 


EDAT1280 

A408 

B 

A409 


ED AT 1 290 


A 

F , 1 72( A ) 

SET UP THt CURRENT LOAD POINT... 

EDATl 300 


LR 

3 , F 

• mm 

EDATl 310 


ST 

3, 32(A) 

• mm 

EDATl 320 


SPACE 

1 


EDATl 330 

A409 

SR 

A, A 

SET UP RA, RB, RE, RF,... 

EDAT 1 340 


LA 

B , S 1 0 0 

■ * • 

EDATl 350 


LA 

£ * A 

• • • 

EDATl 360 


LA 

F , 2 A 

• m • 

EDAT 1370 


LA 

6,1(6) 

CHFCK THE CARD COL. 10 CHAR... 

EDA T 1 3 80 

A410 

CLC 

0( 1,6) , 0 ( B ) 

... 

EDATl 390 


BE 

A500 ( A ) 

BR. IF . FO. AN ALLOWED CHAR. 

EDA T1400 


LA 

B,1(B> 

INCR. RB. 

EDA T 1 410 


RXLE 

A , E , AA 1 0 

LOOP. • • 

EDATl 420 


B 

C110 


EDATl 430 


SPACE 

1 


EDA T 1440 

A500 

B 

BLNK 

BR. OUT AS REO*D... 

EDA T 1450 


B 

CC 

• mm 

EDAT 1 460 


B 

DD 

• mm 

EDAT 1470 


R 

HH 

• mm 

EDAT 1480 


B 

RRR 

• • m 

EDATl 490 


B 

XX 

• mm 

FDAT1500 


TITLE 

•EDATA-C DATA, BLAMK 

AND RETURN PROC. SECTIONS* 

EDAT 1510 


SPACE 

1 


EDATl 520 

* 


* * * 


EDAT 1530 



SPACE 

l 


ED ATI 540 

cc 

01 

Cl02+l*X , F0 i 

ENABLE THG RR*S AT C102... 

E DA T l 550 


01 

C 300+ 1 » X * F 0 * 

« . « C C300 • 

EDA T 1 560 


SPACE 

1 


EDAT1570 

ClOO 

LA 

6,1(6) 

[NCR., SAVE, £ CHECK R6... 

EDAT 1580 


ST 

6, 44(A) 

... 

EO AT l 590 


CR 

6,7 

... 

£ DAT1600 


PH 

BLNK 

BR. OUT IF PAST CARO END. 

EO AT l 6 10 


ClI 

0(6) , X* 70 * 

LOOK FOR THE ' MARK... 

ED ATI 620 


RE 

C300 

...RR. OUT WHEN FOUND. 

EDAT 1630 


SPACE 

1 


EDAT1640 

C 1 02 

RC 

O 

r~ 

o 

c 

** 

EDAT 1650 


SPACE 

l 


ED A T l 660 

C 1 0 3 

C 

9,72(4) 

CHECK FOR PRECEDING ERRORS... 

EDAT 1670 


BE 

ClOO 

...BR. TO AVOID DATA LOADING. 

ED AT1680 


MVC 

01 1,3) ,0(6) 

MOVE A CHARACTER. 

EDA Tl 690 


LA 

3,1(3) 

INCR. R 3 . 

EDAT 1700 


B 

ClOO 


EDAT 1710 


SPACE 

l 


EDA T 1 720 

CUO 

SR 

6,5 

SET UP TROUBLE INOIC... 

EDAT 1 7 30 


LA 

6, 1(6) 

... 

ED ATI 740 


ST 

6,76(4) 

... 

EDAT 1 750 


MV! 

75(4) ,3 

INOIC. TROUBLE. 

EDAT1 760 


B 

RR 


FDAT1770 


SPACE 

1 


EDAT1 780 

C 300 

ec 

15.C400 

** 

EDA Tl 790 


SPACE 

1 


EDAT 1800 


CLC 

0(2,6) ,^X« 7070 * 

CHECK FOR THE ■ 1 MARK. 

EDAT1810 


BNE 

BLNK 

BR. OUT IF NOT A • • MARK. 

EDAT1820 


LA 

6, 1(6) 

INCR. R6. 

EDAT 1 830 


ST 

6,44(4) 

SAVE R6. 

EDAT l 840 


B 

C 1 03 


EDAT 1 850 


SPACE 

1 


EDAT1 860 

C400 

Nl 

C 102+ l , X * OF 1 

DISABLE THE BR. AT C102. 

EDAT 1870 


Nl 

C 300+ l , X 1 OF * 

DISABLE THE BR. AT C300. 

EDAT1880 


B 

ClOO 


EDAT1890 


SPACE 

1 


EDAT 1900 

* 


* * * 


EDAT 1910 


SPACE 

1 


EDAT1920 

BLNK 

MV! 

79(4), 0 

INDIC. NO ERROR. 

EDAT 1 930 


ST 

3,32(4) 

SAVE THE CURRENT LOAD POINT. 

EDAT 1 940 


SPACE 

l 


EDAT1950 

RR 

L 

1,24(4) 

RESET Rl... 

EDAT 1 960 


SPIE 

MF = ( E, (D) 

...TO RESET THE OLD P.I.E. 

ED AT 1970 


L 

0,4(4) 

RESET RD... 

EDAT1980 


LM 

E,C, 12(D) 

...AND RESTORE THE OTHER R*S. 

EDAT l 990 


MV I 

12 ( D ) , X # FF 1 


EDAT 2000 


BCR 

15, e 


EDAT2010 


SPACE 

i 


EOAT202O 

RRR 

MV I 

75(4), 2 

INDIC. RETURN. 

EDAT2030 


R 

RR 


EDAT2040 


TITLE 

• EO AT A-D DATA PROC. 

SECTION* 

EDAT2050 


SPACE 

1 


EDAT2060 

* 


* * * 


EDAT2070 


SPACE 

1 


EDAT2080 



DO 

SR 

A, A 

SET UP R A £ RB . • . 

EDA T2090 


SR 

B,B 

• • • 

EDAT2L00 


XC 

CL 16, CL16 

ZERO HOP, LOP, HOLD, & EXP. 

EDAT2110 


M VC 

TCW,=X'4EDFFFFF' 

SET UP TCW. 

EDAT21 20 


SPACE 

2 


EDAT21 30 

D200 

LA 

6,1(6) 

[NCR. , SAVE, & CHECK R6... 

FDAT2140 


ST 

6, 44(4) 

• • * 

EDAT21 50 


CR 

6,7 

• • • 

EDAT2160 


BH 

08 10 

...SR. OUT IF PAST CARD EMD. 

ED AT 2 170 


SPACE 

I 


E0AT2180 


CLI 

0(6), C'O' 

CHECK FOR NUMERICS... 

EDAT2190 


BL 

osoo 

...RR. OUT IF .LT. 0. 

EDAT2200 


CLI 

0 ( 6 ) ,C * 9' 

... 

EDAT2210 


BH 

Clio 

...BR. TC TROUBLE IF .GT. 9. 

EDAT2220 


MVN 

0(1, 8), 0(6) 

MOVE THE NUMERIC TO 'HOLD'. 

EDAT2230 


LA 

A, 1(A) 

I NCR. THE DIGI T COUNT. 

E DAT2240 


SPACE 

1 


EDAT2250 


TM 

TCI, X'80' 

(0,0) 

EDAT2260 


BO 

D22 0 


EDAT2270 


CH 

A , =H 1 8 1 

CHECK THE DIGIT COUNT... 

EDAT2280 


BL 

D220 

...RR. AROUND IF .LT. 8. 

EDAT2290 


01 

TC 1 , X ' 80* 

DISCONT. THE DIGIT COUNT CHECK. 

E0AT2300 


N I 

TC 3 , X* EF 1 

ALLOW LONG REAL STORAGE. 

EDAT2310 


Ni 

TC 2 , X 1 FD* 

FOR8I0 AN E IN A LONG REAL. 

EDAT2320 


SPACE 

1 


EDAT2330 

0220 

TM 

TC 1 , X 1 40* 

(1,1) 

EDAT2340 


BO. 

0300 


EDAT2350 


LA 

B , 1 ( B ) 

1NCR. THE FRACT. DIGIT COUNT. 

EDAT2360 


SPACE 

2 


EDAT 2370 

D300 

TM 

TC 1 , X * 20 1 

(2,0) 

EDAT2380 


BO 

D400 


EDAT2390 


L 

E , LOP 

BUILD UP THE LONG PRIMITIVE... 

FDAT2400 

D310 

LR 

F,E 

m m m 

EDA T2410 


M 

E , =F 1 10' 

• • • 

EDAT2420 


A 

F.HOLO 

• • • 

EDAT2430 


TM 

TC 1 , X * 10* 

( 3,0) 

EDAT2440 


BO 

0330 


EDAT2450 


ST 

F , LOP 

• • • 

EDAT2460 


N I 

LOP, X' 03' 

• • ♦ 

EDAT2470 


SR 

E,E 

• • • 

EDAT2480 


SLDL 

E , 6 

• • • 

EDAT2490 


ST 

E , HOLD 

• • • 

EDAT2500 


L 

E,HOP 

• • • 

EDAT 2 510 


01 

TCI, X' 10' 

...'OR' BIT 3. 

EDAT2520 


R 

D3 10 


E OAT 2530 


EJECT 



EDAT2540 

0330 

ST 

F , HOP 

• • • 

EDAT2550 


NI 

TClfX'EF* 

...ZERO BIT 3. 

EDAT2560 


R 

0200 


EDA T2570 


SPACE 

2 


EDAT2580 

DA 00 

CH 

A , =H • 2 ' 

CHECK THE DIGIT COUNT... 

EDA T2590 


RH 

C110 

...OR. TO TROUBLE IF .GT. 2. 

EDAT2600 


L 

F , EXP 

BUILD UP THE FXPONFNT • • • 

EDAT2610 


M 

E, =F 1 10' 

... 

EDAT2620 


A 

F , HOLD 

• • • 

EDAT2630 


4 ^ 

ON 



ST 

F , EX P 


B 

D200 


SPACE 

2 

D500 

CLI 

0(6) ,C # ,' 


BE 

D600 


SPACE 

1 


SR 

ltl 


LA 

C, S200 


LA 

E , A 


LA 

F , 20 

D510 

CLC 

0(1, 6), 0(C) 


BE 

D700( 1) 


LA 

C, 1(C) 


BXLE 

1,E,D510 


B 

C110 


SPACE 

2 

0600 

TM 

TC 1 , X 1 08 ' 


BO 

INT 


SPACE 

1 


TM 

TC l , X 1 OA 1 


BO 

D62 0 


SPACE 

1 


A 

B,EXP 

D610 

Nl 

TC3,X'BF* 


SPACE 

1 

0611 

SLL 

B , 3 


ST 

B , EXP 


B 

FLT 


SPACE 

1 

0620 

S 

B , EXP 


R M 

0630 


B 

0610 


SPACE 

1 

0630 

LPR 

B , B 


B 

D611 


SPACE 

2 

0700 

B 

0810 


R 

0820 


B 

0830 


B 

D8A0 


B 

D850 


B 

EJECT 

0860 

D810 

NI 

TC 3 , X* 7F f 


B 

0600 


SPACE 

l 

0820 

TM 

TCI, X* 02* 


PO 

D 8 2 1 


B 

C1I0 

0821 

NI 

TC 1 , X f HO* 


NI 

TC 2t X * 6F • 

0825 

NI 

TC I , X ’ F 7* 


B 

0200 


CHECK FOR A COMMA. . . 
• . . BR. OUT IF COMMA. 


SET UP RC. RE « C RF . . . 


• • • 

CHECK FOR OTHER CHARS... 
...BR. OUT AS REQ* D. 

I NCR. RC. 

LOOP... 

...OR BR. TO TROUBLE. 

A,l) 


5,1) 


COMP. THE SCALING FACTOR... 
...ENABLE DOWNSCALING. 


BLANK. 

+ 

D 

E 

ZERO BIT 16. 


( 6,1 ) 


ALLOW TRACT. DIGIT COUNTING... 
...FORBID 2 .’S IN THE P-PART. 
FORBID + AFTER THE DEC. PT... 
...ALSO, - AFTER THE DEC. PT . 
ALLOW FLOATING. 


EDAT26A0 
EDAT2650 
EDAT2660 
EDA T2 670 
EDAT2680 
EDAT2690 
E DAT 2 7 00 
EDAT2710 
EDAT2720 
EDAT2730 
EDAT27A0 
EDAT2750 
EDAT2760 
EDAT2770 
EDAT2780 
EDAT2790 
EDAT2&00 
EDAT281 0 
EDAT2820 
EDAT 2830 
EDAT28A0 
EDAT2850 
EDAT 2 860 
EDAT2870 
EDAT2880 
EDAT 2 890 
E0AT29D0 
EDAT2910 
EDAT2920 
ED AT 2930 
EDAT29A0 
EDAT2950 
EDAT2960 
EDAT2970 
EDAT29B0 
EDAT2990 
EDAT3000 
EDAT301 0 
EDAT3020 
EDAT3030 
EDAT30A0 
EDAT 3050 
EDAT3060 
EDAT3070 
EDAT3080 
EDAT3090 
EDAT 3100 
EDAT3110 
EDAT3120 
♦EDAT 3130 
EDA T31 AO 
♦EDAT 31 50 
EDAT3160 
EDAT3170 
EDAT3180 


I 


'■mini 




SPACF 

1 


EDAT3190 

D830 

TM 

TCl.X'OP 

17,0) 

EDAT3200 


BO 

D834 


EDAT321 0 


TM 

TC 2 • X 1 BO' 

(8,1) 

E0AT3220 


BO 

0832 


EDAT 32 30 


R 

cno 


EDAT3240 

D832 

N I 

TC? , X 1 6F 1 

FORBID TWO + 'S IN THE P-PART... 

* EDAT3250 




...ALSO, + & - IN THE P-PART. 

EDAT 3260 

0833 

CH 

A , =H* 0* 

CHFCK THE DIGIT COUNT... 

E0AT32 70 


BH 

C110 

...BR. TO TROUBLE IF ,NE • 0. 

EDA T3280 


B 

D200 


E0AT3290 


SPACE 

1 


EDAT3300 

0834 

TM 

TC2 , X 1 40 1 

(9, l ) 

EDAT3310 


BO 

0835 


EDAT 3320 


B 

C110 


FDAT3330 

D835 

N I 

T C2 » X 1 B 7 1 

FORBID TWO +'S IN THE E-PART... 

♦ EDAT 3340 




...ALSO, + & - IN THE E-PART. 

EDAT3350 


B 

0833 


EDAT 3360 


SPACE 

1 


EDAT3370 

0840 

TM 

TC2 » X* 20* 

(10,0) 

FDAT3380 


BO 

D843 


EDAT 3 390 


TM 

TC 2 * X 1 10' 

(11,1) 

EDAT3400 


BO 

D84 2 


ED AT 341 0 


B 

cuo 


F0AT3420 

0842 

N I 

TC2 , X 1 6F * 

FORBID - fi + IN THE P-PART... 

♦E0AT3430 




...ALSO, TWO -• S IN THE P-PART. 

EDAT 3440 


N I 

TC 2 f X 1 FE 1 

ALLOW NEG. INTEGER. 

EDAT3450 


N I 

TC 3 t X 1 OF * 

ALLOW NEG. REAL. 

EDAT 3460 


B 

0833 


EDAT3470 


SPACE 

1 


EDAT3480 

D043 

TM 

TC2 » X 1 08 1 

(12,1) 

EDAT 3490 


BO 

0844 


EDAT3500 


B 

CllO 


E0AT3510 

0844 

N I 

TClfXTB* 

ALLOW NEG. EXP. 

EDAT 3520 


NI 

TC2 , X* B7' 

FORBID - & + IN THE E-PART... 

* EDA T3530 




...ALSO, TWO -*S IN THE E-PART. 

EDAT3540 


B 

D833 


EDAT 3550 


EJECT 



EDAT 3560 

D850 

TM 

TC2 » X * 04* 

(13,1) 

EDAT 3570 


BO 

0851 


EDAT3580 


B 

CllO 


EDAT3590 

0851 

NT 

TC2,X*F9« 

FORBID TWO D'S IN AN ITEM... 

*EDAT3600 




...AND 0 £ E IN AN ITEM. 

EDAT3610 


NI 

TC 3 * X' FF • 

ALLOW LONG REAL STORAGE. 

EDAT 3620 

D852 

01 

TC 1 » X ' 60' 

STOP THE FRACT. DIGIT COUNT... 

♦EDAT3630 




...AND ALLOW EXP . COMP. 

EDAT 3640 


01 

TC l , X * 0 1 • 

SET IJP EXP. SIGN CHECKS... 

EDAT3650 


or 

TC 2 * X* 20* 

• • • 

EDAT 3660 


NI 

TC 1 » X' FD* 

FORBID . IN E-PART. 

EDA T 3670 


SR 

A, A 

ZERO RA. 

EDAT3680 


B 

0825 


EDAT3690 


SPACE 

1 


EDAT3700 

D860 

TM 

TC2 , X * 02* 

( 14 F 1) 

EDAT 3710 


BO 

D861 


EDAT3720 


B 

CllO 


EDAT 3 730 


00 


0861 


IN20 


IN30 


FLT 


FL50 


RL 10 


RL11 


RL 13 


N I 

TC2,X'F9» 

B 

D852 

SPACE 

1 

LM 

E , F , HUP 

CH 

E , = H 1 31' 

BH 

C110 

SLL 

F , 6 

SRDL 

E,6 

TM 

TC 2 , X* 01* 

BO 

I N20 

LNR 

F,F 

C 

9,72(4) 

BE 

I N 30 

ST 

F , 0 ( 3 ) 

LA 

3,4(3) 

SPACE 

1 

TM 

T C 3 , X 1 80* 

BO 

DD 

B 

BLNK 

SPACE 

1 

LD 

2,=D*0.0» 

LM 

A , B , HOP 

LTR 

E , A 

BZ 

FL 50 

SR 

F , F 

SRDL 

E, 10 

STM 

E,F, HUP 

01 

HOP , X 1 4F * 

AD 

2 , D WD 

LTR 

B,B 

BZ 

RL 10 

SR 

A, A 

STM 

A , B , HUP 

01 

HOP » X # 4E 1 

AD 

EJECT 

2 , D WD 

SPACE 

2 

LTDR 

2,2 

RZ 

RL 30 

L 

C , EXP 

LTR 

C,C 

BZ 

RL20 

SPACE 

1 

L 

1 , =A l LFH3 

TM 

TC 3, X • 40 1 

BO 

RL 15 

SPACE 

l 

CH 

C * -H f 608 • 

BL 

RL13 

SH 

C , 600* 

DD 

2,600(1) 

B 

RL 11 

DO 

2 , 0 ( C , 1 ) 

R 

RL20 


FORBID TWO F 1 S IM AN ITEM... *E0AT3740 
...AND E & D IN AN ITEM. EDAT3750 

EO AT 3 760 
EDAT3770 

LOAD THE LONG PRIMITIVE... EDAT3730 

...£ CHECK THE H-ORDER PART. EDAT3790 

...BR. TO TR0U8LE IF TOO LARGE. EDAT3800 

LINK UP THE PARTS... EDAT3810 

...£ R-SHIFT IT ALL INTO RF. EDAT3820 

( 15 » l ) E0AT3830 

EDAT3840 

MAKE THE RESULT NEG. , IF REQ 1 D. EDAT3850 

CHECK FOR PRECEDING ERRORS... EDAT3860 

...BR. TO AVOID DATA LOADING. EDAT3870 

STASH THE INTEGER WHERE R E 0 1 D . EDAT3880 

INCR. R3. EDAT 3890 

EDAT3900 

(16,1) EDAT 3910 

EDAT3920 
FDAT3930 
EDAT 3940 

FLOAT THE LONG PRIMITIVE... EDAT3950 

... EDA T3960 

EDAT3970 
EDA T3980 

... EDAT3990 

... EDAT4000 

... EDA T40 1 0 

... EDAT4020 

... EDAT4030 

... E0AT4040 

... EDAT 4050 

... EDAT4060 

... EDAT4070 

... EDAT4380 

... EDAT4090 

EDAT4100 
EDAT 4110 

CHECK THE FLOATED RESULT... EDAT4120 

...BR. OUT TO AVOID SCALING. EDAT4130 

...£ SCALE AS REQ v D. . . EDAT4140 

EDAT4150 

EDAT4160 

EDAT4170 

... EDAT4180 

(17,1) EDAT4190 

... EDAT4200 

EDAT4210 

... EDAT 4220 

... EDAT4230 

EDAT A 240 
EDAT 4250 
EDAT4260 

... EDAT4270 

EDA T4280 




SPACE 

1 


EDAT4290 

RL15 

CH 

C,=H'600' 

• . • 

EDAT4300 


BH 

CLIO 

•••OR. TO TROUBLE IF TCQ MUCH. 

EDAT4310 


MD 

2,0<C,1) 

... 

EDAT4320 


SPACE 

1 


EDAT4330 

RL20 

TM 

TC3,X'20» 

(13,1) 

EDAT 4340 


R0 

RL 30 


ED AT 4350 


LNDR 

2 ? 2 

MAKE IT MEG. IF REQ'D. 

EDAT4360 


SPACE 

1 


E0AT4370 

RL30 

TM 

TC 3 ? X 1 10* 

(19,1) 

EDA T4380 


BO 

RL50 


EDAT 4 390 

RL 31 

C 

9, 72(A) 

CHECK FOR PRECFDING ERRORS... 

EDAT4400 


RE 

mo 

• ••HR. TO A VO I 0 OATA LOADING. 

EDAT 4410 


STD 

2,0(3) 

STASH THE L-REAL WHERE REG'D. 

EDAT 4420 


LA 

3, 8{ 3) 

INCR. R3. 

EDAT4430 


R 

I N 30 


EDAT4440 


SPACE 

2 


EDAT4450 

RL 50 

LTDR 

2,2 

RUUND OUT AS RED* 0* • • 

EDAT4460 


BZ 

RL 5 5 

• • • 

EDAT4470 


STD 

2 , DWD 

• * • 

EDAT4480 


NC 

HOP , =X ' FFO 00000* 

• • • 

EDAT 4490 


NC 

LOP , =X* BOGOOCOO ' 

• • • 

EDAT4500 


AD 

2, DWD 

• • • 

EDAT451 0 

RL 55 

C 

9,72(4) 

CHTCK FOR PRECEDING ERRORS... 

EDAT4520 


RF 

If4 30 

...RR. TO AVOID DATA LOADING. 

EDAT 4530 


STE 

2,0(3) 

STORE THE S-REAL WHERE REO'D. 

E0AT4540 


LA 

3, 4 ( 3 ) 

INCR. R 3. 

EDAT4550 


B 

I N 30 


EDAT4560 


TITLE 

'EDATA-H DATA PROC 

. SECTIUN' 

EDAT4570 

* 


* * * 


EDAT4580 


SPACE 

1 


FDAT4590 

HH 

SR 

A, A 

SET UP R A • • • 

FDAT4600 


XC 

CL 8 , CL 8 

ZERO LOP & HOLD. 

EDAT 461 0 


01 

H510+i,X'r0' 

FNABLE THE BR. AT H510. 

EDAT4620 


SPACE 

I 


EDAT 4630 

H200 

LA 

6,1(6) 

INCR* , SAVE, & CHECK R6... 

EDAT 4640 


ST 

6,44(4) 

• • • 

EDAT 4650 


CR 

6,7 

• • » 

E0AT4660 


BH 

H600 

...BR. OUT IF PAST CARD END. 

EDAT 46 70 


SPACE 

1 


EDAT4680 


CLI 

0(6), C',* 

CHFCK THE DATA FIELD CHAR... 

EDAT4690 


BL 

H 300 

...RR. GUT IF .LT. 'COMMA' . 

EDAT4700 


BE 

H500 

...DR. OUT IF .EQ. 'COMMA' . 

EDAT4710 


CLI 

0( 6), C' O' 


EDAT4720 


PL 

CllO 

...RR. TO TROUBLE IF .LT. 0. 

EDAT4730 


CLI 

0(6) ,C' 9' 


EDAT 4740 


RH 

CllO 

...BR. TC TROUBLE IF .GT. 9. 

E DAT4750 


MV N 

0( 1,8), 0(6) 

MOVE THF NUMERIC TO 'HOLD'. 

EDAT4760 


LA 

A, 1(A) 

INCR. RA. 

EDAT4770 


SPACE 

1 


EDAT4780 


L 

F , LOP 

BUILD UP THE H- FORMATTED DATA... 

EDAT 4 790 


M 

E , = F • 10* 

... 

EDAT4800 


A 

F, HOLD 

. • • 

EDAT4810 


CH 

F , = H' 32767' 

... 

EDA T4820 


RH 

CllO 

...RR. TO TRGUDLE IF .GT. 32767. 

EDAT4Q30 



Lh 

O 



ST 

F f LQP 

• • « 

EDAT 4840 


B 

H200 


EDAT 4850 


SPACE 

l 


EDAT4860 

H300 

CLl 

016) ,C ' * 

RECHECK THE DATA FIELD CHAR... 

EDA T4870 


BE 

H600 

...BR. OUT IF BLANK. 

EDAT4880 


CU 

0 ( 6 ) ,C * +* 

• • • 

EDAT4890 


BE 

H410 

...OR. OUT IF PLUS SIGN. 

EDAT4900 


CLl 

0 ( 6 } * C * -* 

• • ♦ 

EDAT4910 


BE 

H400 

...BR. OUT IF MINUS SIGN. 

EDAT 4920 


B 

C110 

...BR. TO TROUBLE OTHERWISE. 

EDAT4930 


SPACE 

1 


EDAT 49^0 

H400 

Nl 

H5 10+ l » X' OF 1 

DISABLE THE BR. AT H510. 

EDAT4950 

H410 

CH 

A , ~H* 0* 

CHECK FOR SIGN EMBEDMENT... 

EDAT 4960 


BH 

C110 

...BR. TO TROUBLE IF EMBEDDED. 

EDAT4970 


B 

H200 


EDAT4 980 


SPACE 

1 


EDAT4990 

H500 

C 

9,72(4) 

CHECK FOR PRECEDING ERRORS... 

EDAT5000 


BE 

H530 

...BR. TO AVOID DATA LOADING. 

EDAT50I0 


LH 

B,L0P+2 


EDAT 5020 

H510 

BC 

1 5 t H52 0 

** 

EDA T5030 


LNR 

B , B 

SET SIGN MINUS IF REQ*D. 

EDAT5040 

H520 

STH 

B , 0 ( 3) 

STORE THE H- INTEGER. 

EDAT5050 


LA 

3,2(3) 

INCR. R3. 

EDAT5060 

H530 

BC 

15, HH 

** 

EDAT5D70 


ni 

H530 + l , X* FO * 

ENABLE THE BR. AT H530. 

EDAT5080 


B 

BLNK 


E0AT5090 


SPACE 

1 


EDAT5100 

H600 

Nl 

H530+ 1 , X * OF • 

DISABLE THE BR. AT H530. 

EDAT5110 


B 

H500 


EDAT5120 


TITLE 

•EOATA-X DATA PROC. 

SECTION* 

EDAT5130 


SPACE 

1 


EDAT5140 

♦ 


* * * 


EDAT5L50 


SPACE 

1 


EDAT5I60 

XX 

SR 

A, A 

ZERO RA... 

SDAT5170 


xc 

CL8 , CL8 

ZFRO LOP £ HOLD. 

EDAT5180 


SPACE 

1 


EDAT 5190 

X200 

LA 

6,1(6) 

INCR., SAVE, £ CHECK R6... 

EDAT5200 


ST 

6,44(4) 

• • • 

EDAT5210 


CR 

6,7 


EDAT5220 


BH 

X500 

...BR. OUT IF PAST CARD END. 

EDAT5230 


SPACE 

l 


EDAT5240 


CLl 

0(6), C*, 1 

CHECK THE OATA FIELD CHAR... 

EDAT5250 


BE 

X400 

...BR. OUT IF COMMA. 

EDAT5260 


CU 

0(6), CO 1 

• t « 

EDAT5270 


BL 

X 300 

...BR. OUT IF .LT. 0. 

EDAT5280 


CLl 

0 ( 6 ) , C • 9* 

... 

EDAT5290 


BH 

Clio 

...BR. TO TROUBLE IF .GT. 9. 

E0AT5300 


MVN 

0( 1,8), 0(6) 

MOVE THE NUMERIC TO •HOLD*. 

EDAT5310 


L 

B, HOLD 

SET UP RB. 

EDAT5320 


SPACE 

1 


EDAT5330 

X220 

LA 

A, It A) 

INCR. £ CHECK R8... 

EDAT5340 


CH 

A , =H * 8* 

• • • 

EDAT5350 


BH 

CUO 

...BR. TO TROUBLE IF .GT. 8. 

E0AT5360 


SPACE 

i 


EDAT5370 


L 

C , LOP 

BUILD UP THF X. FORMATTED DATA... 

EDAT5380 




SLL 

C,4 

• • • 




EDA T5390 


OR 

CtB 

• • • 




EDAT5400 


ST 

C t LOP 

• mm 




EDAT5410 


R 

X200 

m • • 




EDAT5420 


SPACE 

1 





EDAT5430 

X300 

CLI 

0 ( 6 ) * C 1 1 

RFCHECK THE DATA 

FIELD CHAR... 

EDAT5440 


RE 

X500 

...BR. 0U1 

IF BLANK. 


EDAT5450 


CLI 

0 ( 6 ) » C • A* 

• mm 




EDAT5460 


BL 

cno 

. • • B R • TO 

TROUBLE IF 

.LT. • A* . 

EDAT5470 


CLI 

0(6) f C , F' 





EDAT5480 


BH 

cuo 

• • • 8R. TC 

TROUBLE IF 

.GT. 1 F* • 

EDAT5490 


MVN 

0( 1 1 8 ) * 0 ( 6 ) 

MOVE THE NUMERIC 

TO • 

HOLD*. 

EDAT5500 


L 

B, HOLD 

PICK UP THF NUMERIC.. 

. 

EDAT5510 


AH 

B , =H * 9 1 

...6 REMOVE ITS 

BIAS. 


EDAT5520 


R 

X220 





EDAT5530 


SPACE 

1 





EDAT5540 

X400 

C 

9,72(4) 

CHECK FOR 

PRECEDING ERRORS... 

EDAT 5550 


RE 

X410 

...RR. TC 

AVOID 

DATA 

LOADING. 

EDAT5560 


SPACE 

1 





EDAT5570 


L 

B, LOP 

STASH THE 

DATA WHERE 

REQ'D... 

EDAT5580 


ST 

B, 0(3) 

■ • ♦ 




EDAT5590 


LA 

3,413) 

INCR. R3. 




EDA T5600 


SPACE 

1 





EDAT561 0 

X410 

BC 

15, XX 

** 




FDAT5620 


01 

X410+1,X , F0* 

ENABLE THE 

BR. AT X410. 

EDA T5630 


B 

BLNK 





E0AT5640 


SPACE 

1 





EDAT5650 

X500 

N 1 

X4 1 0+ 1 , X * OF • 

DISABLE THF BR. 

AT X410. 

ED AT 5660 


B 

X400 





EDAT5670 


TITLE 

1 EOAT A-ERAS I BLE 

STORAGE AND CONSTANTS. 1 



EDAT5680 


SPACE 

1 





EDAT5690 

DWD 

DS 

OD 





FDAT5700 

CL 16 

DS 

0CL16 





EDAT5710 

HOP 

DS 

F 





E0AT5720 

CL8 

DS 

OCL 8 





EDAT5730 

LOP 

DS 

F 





EDAT 5740 

HOLD 

DS 

OF 





FDAT5750 


DS 

CL3 





EDAT5760 

LOB 

DS 

C 





EDAT5770 

EXP 

DS 

F 





EDAT 5780 


SPACE 

1 





EDAT5790 

TCW 

DS 

OF 





FDAT5800 

TCI 

DS 

C 





FDAT58 1 0 

TC2 

DS 

c 





EDAT5820 

TC3 

DS 

C 





EDAT5830 

TC4 

DS 

c 





EDAT5840 


SPACE 

1 





EDAT5850 

SlOO 

DC 

c* • 





FDAT5860 


DC 

C'C' 





EDAT5870 


DC 

C ' 0 1 





EDAT5880 


DC 

C 1 H 1 





EDAT5890 


DC 

C * R * 





EDAT5900 


DC 

C* X* 





EDAT 591 0 


SPACE 

1 





EDAT5920 

S200 

DC 

C # ■ 





EDAT5930 



DC C* 

% • 

EDAT5940 


DC C' 


EDA T5950 


DC C* 

_ • 

EDAT5960 


DC C 11 

IV 

EDAT5970 


DC C f 

E* 

EDAT5980 


SPACE I 


EDAT5990 


LTORG 


EDAT6000 


TITLE ' EOiTA-CSFCTS LFH2 AND LFH3. ' 

EQAT6Q10 


SPACE 1 


EDAT6020 

* 

SPACF L 

********* 

FDAT6030 
EDAT 6040 

LFH2 

C SEC T 


EDAT6050 

FWO 

DS F 

+0 

EDAT6060 

RD 

OS F 

+ 4 

EDAT 6070 

RZ 

DS 4F 

+ 8 

EDAT6080 

R1 

OS F 

+24 

EDAT6090 

R2 

DS F 

+ 28 

EDAT 6100 

R3 

DS F 

+ 32 

EDAT6110 

R4 

DS F 

+ 36 

EDAT 6120 

R5 

DS F 

+ 40 

EDAT 6130 

R6 

DS F 

+44 

EDAT 6140 

R7 

DS F 

+ 48 

EDAT 6150 

R8 

DS F 

+ 52 

EDAT 6160 

R 9 

DS F 

+ 56 

E0AT6170 

RY 

DS 3F 

+60 

EDAT6180 

N 1 

DS F 

+ 72 

EDAT 6190 

N 2 

DS F 

+ 76 

E0AT6200 

N3 

DS F 

+ 80 

EDAT6210 

N4 

DS F 

+ 84 

EDAT6220 

CARD 

DS CL 80 +88 

E0AT6230 

BLANK 

DC CL4 1 1 

E0AT6240 

R30 

DS F 

+ 172 

EDAT6250 

PSW 

DS 2F 

+ 176 

EDAT6260 

K3 

DS 2H 

+ 184 

EDAT6270 

CBLOCK 

DS 4400F +188 

EDA T 6280 


SPACE t 


EDA T6290 

* 

SPACE L 

********* 

F0AT6300 

CDAT6310 

LFH3 

C SECT 


EDAT6320 

* 


TABLE OF 10**N. 

EDAT 6 330 

* 



EDAT 6 340 

TENS 

DC D* 

l.OE+O, l.OE+L, l.OE+2, 1.0E+3, l.OF+4, l .OE+5* 

EDAT6350 


DC D* 

1.0t+6,l.0E+7,1.0E+8,l.0E+9,l.0E+10' 

EDAT6360 


DC D* 

l. 0E+ 1 1, l. CE + 12, 1.0E+ 13,1.0 6 + 14,1 .OE+15* 

EDAT6370 


DC IV 

l.0E+lA,1.0E+t7,l.0E+13,l.0£+l9,l.0E+20' 

EDAT6380 


EJECT 


EDAT6390 


DC D* 

L.OE+21, l.OF+22, 1.0E+23, 1.0F+24, I.OE+25’ 

EDAT6400 


DC D' 

1.0E+26 , l.OE+27,1. OE + 2 8,1.0 E+29,1 .06+30' 

EDAT6410 


DC D' 

L.0E + 3L , l.OE+32,1.0E+33,1.0E+34,l .OE+35' 

EDAT 6420 


DC O' 

1.0E+ 36, 1.0E+ 37, l.OE+38, 1.06+39,1.06+40* 

EDAT 6430 


DC O' 

1.0E+4l,l.OE+42,1.0i£+43,1.0E+44,l .OE+45* 

EDAT6440 


OC O' 

l.OE+46, l.OE+47 ,1.05+48, l.OE+49,1 .06+50* 

EDAT 6450 


DC D * 

1.0E+5l,1.0E+52,1.0E+53,l.0E+54,l.0£+55' 

EDAT 6460 


DC D* 

l.OE+56, l.OE +57,1. OE+58, 1.0 E+59, l.Ot+60' 

EDAT6470 


DC D* 

l.OE +61 »l.CE+62,1.0E+63, l.OE +64, 1.06+65* 

EDAT6480 



■tt- ir> t'l 


OC D» L.OE+66, l.OE+67, 1.0E + 68, L. OE+69,1. 0E+70' 

DC D* l.OE+71, 1.0E+72,1.0E+73,1.0E+74, l.Ot+75* 

FND 

/* 

// EXEC ASSEMBLR^ARM^LOAD.DECK' 

/ / SQURCF . SYS I N OD * 

SPIE TITLE '...DIAGNOSTIC SUPPRESSOR SUBPROGRAM...' 

* 9/17/66 - L.F.H. 

* MOD. 3/07/69 - L.F.H. 

4 

* USAGE. ..CALL SPIE( LOW ) 

* WHERE LDW IS A PARAMETER CALLING FUR SETTING A NEW P.I.E. OR RE- 

* SETTING THE OLD P.I.E. DEPENDING ON WHETHER IT IS .GT. ZERO OR .LT. 

* ONE, RESPECTIVELY. IF THE NEW P.I.E. HAS BEEN SET, IT REMAINS IN 

* EFFECT ON SUBSEQUENT CALLS ON SPIE WITH LDW .GT. ZERO*, THE OLD P.I.E. 

* IS RESET IF A SUBSEQUENT CALL ON SPIE HAS LDW .LT. ONE. 

* ...GOOD LUCK-L.F.H. 

* 

SPIE CSECT 

ENTRY YECH , R 1 

USING * , l 5 

R S100 

DC CL6' 5SPIE ' 

AREA DS 1 OF 

100 STM 14.L2.12U3) 

LA 12, SAR E A 

ST 13,4(12) 

ST 12,6(13) 

LR 13,12 

4 

BALR 2,0 

USING * » 2 

DROP 15 

4 

L 3,0(1) 

L 4,0(3) 

C 4 , =F 1 1 • 

BL S 300 

4 

L 4, FLAG 

C 4, =F * 1 * 

RE S200 

MVC FLAG,=F’l» 

4 

SPIF YECH, ((1,15)) 

ST l . R L 

4 

S200 L 13,4(12) 

LM 14,12,12(13) 

MV I 12(13), X'FF' 

YECH BCR 15,14 

4 

S300 L 4, FLAG 

C 4 , =F 1 0 * 


EDAT6490 
EDAT 6500 
EDAT 6510 


SPIF0010 
SPIE0020 
SP I EGO 30 
SPIE0040 
SPIE0050 
SPIE0060 
SPIF0070 
SPIE0080 
SP I F 0090 
SPIEOLDO 
SPIEOUO 
SP I EO 120 
SPIE0L30 
SPIE0140 
SP I EO 150 
SP I FO 160 
S PI E01 70 
SPIE0180 
SP I EO 190 
SPIE0200 
SPIE0210 
SPIE0220 
SP T F0230 
SPIE0240 
SPIE0250 
SPIE0260 
SP I E 0270 
SPIE0280 
S P I EO 290 
SP I EO 300 
SP I EO 310 
SP I EO 320 
SPIE0330 
SPIE0340 
SPIFO350 
SPIE0360 
SPIE0370 
SPIE03B0 
SPIE0390 
SPIE0400 
SPIF0410 
SPIE0420 
SPIE0430 
SPIF0440 
SPIE0450 
SPIF0460 
SPIE0470 
SPIE04B0 
SPIE 0490 


Ln 

CO 




RE S 200 

SP I EO 500 


L 1? R1 

SPIF0510 


SPIE MF= ( E * (1)) 

SPIE0520 


XC FLAG, FLAG 

SPIF0530 


B S 200 

SPIE0540 

* 


SPIF0550 

R1 

DS F 

SPIE0560 

FLAG 

DC F ' 0 ' 

SPIE0570 


LTORG 

SPIE0580 

/* 

END 

SPIE0590 

// EXFX 

ASSEMBLR* PARM= • LOAD, DECK ' 


//SOURCE. SYSIM DO * 


STIM 

TITLE '...TIME INTERVAL MEASURING SUBPROGRAM...' 

STIME010 

* 6/14/68 - L.F.H. 

STIME 020 

* MOD. 

8/30/68 - L.F.H. 

STIMF030 

* 


ST I ME 040 

STIME 

C SEC T 

ST IME050 . 


ENTRY TTIME f Tl 

STIME060 


USING *,15 

STIME070 


B SIOO 

STIME080 


DC X ' 05 1 

STIME090 


DC CL5 1 STIME* 

ST IMF 1 00 

SiOO 

STM 14,12,12(13) 

STIMEl 1 0 


LA 15, TT IME-STIMF ( 15) 

STIME120 


USING TT I MF , 15 

ST I ME 1 30 


N I T 1 02+ 1 , X' OF' 

STIME140 


M VC T2 , Ti 

ST1ME150 


B T100 

S T I ME 160 

* 


ST I ME 170 

TTIME 

STM 14, 12, 12(13) 

ST I ME 1 80 


L 3,0(1) 

STIME190 


L 4 , T 1 

STIME200 


01 T102+1 ,X* FO' 

STIME210 

T 100 

LR 2,15 

STIME220 


USING TTIME, 2 

STIME230 


DROP 15 

STIME240 


LA 12.SAKEA 

SIIME250 


ST 13,4(12) 

STIME260 


ST 12,8(13) 

STIME270 


LR 13,12 

STI ME 280 

T102 

DC 0, T200 ** 

5TIME290 

+ 


STIME300 

T 103 

STIMER TASK,RTN,TUINTVL=T2 

ST I ME 31 0 

T 1 05 

L 13,4(12) 

STIME 320 


LM 14,12,12(13) 

ST I ME 330 


BCR 15,14 

ST1ME340 

* 


ST I ME 350 

T200 

TTIMFR CANCEL 

STIME360 


L 4, Tl 

STIME370 


SR 4,0 

ST 1 ME 380 


ST 0 , T 1 

STIME390 


ST 4,0(3) 

ST1ME400 


B Tl 05 

STIME410 

* 


STIME420 



RTN 

ABEND 

4095, DUMP* STEP 


STIME430 


F JECT 



STIME440 

* 




STIME450 

T 1 

DC 

F' 138461538' 

T-UN l T S=60 MINUTES. 

STIME460 

T 2 

DS 

F 


STIME470 

SARFA 

DS 

18F 


STIME480 

/* 

// EXEC 

END 



STIMF490 

ASSEMBLR,PARM=' LOAD, DECK' 



//SUURCF 

- SYSIN 

DD * 



DATA 

TITLE 

•...INITIAL DATA IN CSECTS DATA AND CONS...* 

DAT A0010 

* 




DATA0020 

* 2/21/69 - L 

. • F • H. 


DAT A0030 

* MOO. 

9/03/69 - L.F.H. 


DAT A0Q40 

* 




DAT A0050 

DATA 

C SECT 



DATA0060 

* 




DATA0070 

Z 

DC 

2400* 0.0* *800* ' 


DATA0080 

* 




DATA0090 

* 


...PROGRAM CONTROL PARAMETERS... 

DATA0100 

* 




DATAOUO 

NEO 1 

DC 

F • 3 * NO. OF L ST 

PHASE DERIV. EO'S. 

DAT A01 20 

NEQ2 

DC 

F 1 6 1 NO. OF 2ND 

PHASE DERIV. EQ'S. 

DAT AOL 30 

Jl 

DC 

F • 0* BR. PARM. 

SET BY SURP. SETH FOR SUBP. RK . 

DAT A01 40 

J 2 

DC 

F 1 5 1 1ST PHASE 

PRINT FREG. 

D AT AOl 50 

J3 

DC 

F • 5 * H ADJUST. 

PARM. 

DATA01 60 

J4 

DC 

F • 0 * CASE NO. 


DAT AOl 70 

J 5 

DC 

F ' 0 1 2ND PHASE 

EULER1AN POLE INDICATOR. 

DAT AOl 80 

* 




DAT AOl 90 

* 


...CONVERSION CONSTANTS... 

DAT A0200 

* 




DAT A02 10 

CONS 

C SEC T 



DAT A0220 

TWP I 

DC 

D*6. 28 3 1853071795864769* 

2 *P I 

DAT A0230 

CRTD 

DC 

D*57. 29577951308232' 

RADIANS TO DEGREES 

DATA0240 

CDTR 

DC 

D' 1. 74532 92 5 1 9943 3 OE -2 1 

DEGREES TC RADIANS 

DAT A0250 

PID2 

DC 

D' L. 570796326794897' 

PI/2.0 

DAT A0260 

PID4 

DC 

D' .7853981633974483' 

PI/4.0 

DAT A0270 


END 



DATA0280 

/* 






// EXEC LlNKGO 

//GQ.SYSUDUMP DD SYSOUT=A, SPACF= ( TRK, ( 8)) 
//GQ.DATA5 DD * 



...JAVELIN ROC STUDY DATA 

- PART 

1 . 1 ... 

JROCOOiO 

JR0C0020 


NOMINAL DATA WITH RELEASE 

AT VARYING GAMMA ANGLES. 

JR0C0030 

JR0C0040 

1920 

C*... JAVELIN ROC STUDY PART 

1.1 - 

N/J ALIGNMENT...' 

JR0C0050 

8 

D4. 88281250-4, 9. 500 


HO, 01 . 

JR0C0060 

80 

D. 804167 DO* .4286D0* 1 .45600 


Dl ».)? ,03. 

JRQC0070 

136 

D3.50333D0 


D 8. 

JROC0080 

160 

D.4D0 * . 1 17800, . 3466D0 * . 42D0 


M, JCYl, JCY2, JCY3. 

JRUC0090 

192 

DODO, .032 1900 * ODO 


JCY4, JCY5, JCY6. 

JROCOIOO 

240 

D7.5D0 


JVXL. 

JR0C01 1 0 

328 

D6D 1 * 5D-6 


YT , YF PS . 

JR0C0120 

352 

D515D0 


AX 1 . 

JRDC0130 



NASA-Langley, 1971 


2012 

D4 


J2-PRI NT FREC. 

JR0C0140 

2016 

n 15 


J3-H ADJ. PARM D 

JR0C0150 

177 

R 



J ROC 0 1 60 

2012 

032768 


J2-PR t NT FRECc 

JR0C0L70 

328 

D3D 1 


YT. 

JR0C01 80 

77 

R 



JROC0190 

328 

Dl. 5DL 


YT. 

JROC0200 

77 

R 



JRCC0210 

328 

012. 5D0 


YT. 

JROC022O 

77 

R 



JR OC 0230 

328 

D ID l 


YT. 

JR0C0240 

77 

R 



JRQC0250 





JR0C0260 





JRUC0270 


EFFECT OF CLAMSHELL DENSITY CHANGE 

TO 

l.l OF NOMINAL. 

JRDC0280 





JR0C0290 

160 

D.44D0, . 12 9600, . 38 1300, .462 DO 


M t JCY1, JCY2, JCY3. 

JRDC0300 

192 

DODO, .0354100,000 


JCY4, JCY5, JCY6. 

JROC0310 

328 

D6D1, 5D-6 


YT , YEPS. 

JR0C0320 

2012 

04 


J2-PR I NT FREQ. 

JR0C0330 

3170 

R 



JR0C0340 

2012 

D32768 


J2-PRINT FREQ. 

JR0C0350 

328 

D3D 1 


YT. 

JRQC0360 

70 

R 



JR0C03 7 0 

328 

D1.5D1 


YT. 

JROC0380 

70 

R 



JRQC0390 

328 

012. 5D0 


YT. 

JR0C0430 

70 

R 



JR OC 04 1 0 

328 

D ID 1 


YT. 

JROC0420 

70 

R 



JR0C0430 





JROC044O 


EFFECT OF CLAMSHELL DENSITY CHANGE 

TO 

0.9 OF NOMINAL. 

JR0C0450 





JRUC0460 

160 

D. 3600 , . 106000, . 3 1 1900, . 37 800 


M, JCY1, JCY2, JCY3. 

JR0C0470 

L92 

DODO, .0289700,000 


JCY4 , JCY5 , JCY6 • 

J R0C0480 

328 

D6D 1 


YT. 

JR0C0490 

2012 

D4 


J2-PRINT FREQ. 

JROC0500 

4170 

R 



JROC05 1 0 

2012 

D 3276 8 


J2-PRINT FREQ. 

JR0C0520 

328 

D3D1 


YT. 

JR0C0530 

70 

R 



JR0C0540 

328 

01.5D1 


YT. 

J R0C055O 

70 

R 



JRQCO 560 

328 

D12.5D0 


YT. 

JR0C0570 

70 

R 



JR0C0580 

328 

D ID 1 


YT. 

JROC0590 

00 

R ...END OF JOB INPUT... 



JR OC 0600 

/* 
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